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INTRODUÇÃO AOS CÉREBROS ELECTRÓNICOS 


O computador ABCD 


OS COMPUTADORES. O computador ABCD 
(que não existe mas podia existir) é um compu- 
tador electrónico, de programação interna, aritmé- 
tico e universal, imaginado para simplificar esta 
exposição. Naquela designação acrescenta-se ao 
substantivo computador — que em si quer dizer 
uma grande máquina de cálculo automático — 
vários atributos com que se distinguem os prin- 
cipais computadores existentes de outras máquinas 
de cálculo automático. Do MARK I (o calculador 
automático terminado em 1944 na Universidade 
de Harvard e o primeiro dos grandes computa- 
dores) às máquinas que agora equipam os centros 
de cálculo automático do mundo há já uma longa 
evolução ; aqueles atributos são características do 
estado actual dessa evolução. 

O adjectivo electrónico opõe-se a mecânico e 
electromecânico. Quer dizer que os dispositivos 
com que agora se fabricam us computadores são 
dispositivos electrónicos. Electrónico não significa 
bizarro, caro, frágil ou inseguro; no estado ac- 
tual da tecnologia significa directo, barato, resi- 
liente e seguro. 

O computador de programação interna é um compu- 
tador que não recebe as suas ordens, passo a 
passo, do exterior. O programa que em dado 
cálculo tem de cumprir é armazenado nos seus 
registos no princípio do cálculo, como os pró- 
prios dados sobre que actuará; as ordens são 


ror GUSTAVO DE CASTRO 


Matemático (L.N. E. €,) 


exploradas logo antes de serem cumpridas, no 
decurso do cálculo. 

O computador aritmético é o computador que, 
como as máquinas de secretária, executa as ope- 
rações aritméticas simples, que são a base das 
suas possibilidades. Distingue-se assim dos cal- 
culadores analógicos (de que a régua de cálculo 
e o contador da electricidade são exemplos) por- 
que nos computadores analógicos se não executam 
operações aritméticas, por processos como os que 
as pessoas usam, mas se fazem medições de 
grandezas. Os computadores aritméticos também 
se chamam digitais; os computadores analógicos 
também se poderiam chamar geométricos. A ponte 
de Wheatstone, permitindo, com uma força elec- 
tromotriz e algumas resistências, a realização de 
multiplicações e divisões, é um computador ana- 
lógico. 

Um computador universal opõe-se às máquinas 
construídas para a execução de cálculos de uma 
categoria bem definida, ao serviço de um utiliza- 
dor que os tem de realizar, com novos dados, 
muitas vezes. Tais como as máquinas de secre- 
tária, os computadores universais pretendam 
então ser instrumentos de cálculo matemático, 
aptos a satisfazerem todas as solicitações que 
assumam a forma de um cálculo. Pelo contrário, 
pode querer-se a realização contínua dum mesmo 
programa, sem exigências de precisão mas com 
os resultados numa forma expressiva; um cal- 
culador, não universal mas adequado ao programa 
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imutável, que dê o resultado com uma agulha ou 
num gráfico pode ser preferível. 


A ANÁLISE NUMÉRICA. A afirmação de que 
os principais cérebros electrónicos assentam as 
suas possibilidades nas quatro operações deve 
completar-se dizendo que as fazem a velocidades 
muito grandes, com pequena incidência de erros, 
raros e curtos tempos de avaria. Para as pessoas 
sem grande experiência pode parecer estranho 
que, numa ciência em grande parte ocupada na 
determinação de números, tudo se reduza às 
quatro operações. Esclareçamo-nos com um exem- 
plo. 

O estudo dos integrais definidos é um dos 
problemas mais importantes da análise mate- 
mática; o seu cálculo, nos termos da análise 
infinitésimal, é um problema muito difícil que o 
principiante considera razoavelmente o mais for- 
midável dos que se lhe podem propor. O problema 
é, nos termos em que se põe, quase sempre 
insolúvel se não houver cuidados na sua formu- 
lação (circunstância que escapa, aliás. ao princi- 
piante). Como problema de análise numérica, 
porém, o cálculo de um integral definido não é um 
problema difícil. 

A análise numérica é precisamente a disciplina 
onde se criam e arrumam as técnicas pelas quais 
os cálculos matemáticos se reduzem às quatro 
operações. Seja, por exemplo, o cálculo de um 
integral definido. 


fia fi 6) ds (1) 


em que a e b são números e f(x) a função a 
integrar, representada no intervalo (a,b) como 
se vê na Figura 1. 


8 
Y Af 


0 q q; b X 


Fig. 1 — Regra do trapézio 


Segundo a teoria, 1 é a área limitada pela 
curva, o eixo e as ordenadas levantadas em a 
e b. Os métodos numéricos são em geral méto- 
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dos aproximados; o problema que se põe é o da 
escolha da superfície, cuja área seja fácil de cal- 
cular e pouco diferente da área I. Uma solução 
simples está em tomar a área do trapézio, abBA ; 
teremos então 


[ima=i t+) +E (2) 


e o que se sugere é que se despreze E. O erro 
que se comete é o valor de E, isto é, a área limi- 
tada pela curva e a corda AB, é um erro consi- 
derável. 

Sem abandonar o princípio, podemos agora 
levantar uma ordenada por a,, a meia distância 
entre a e b, e aplicar às duas partes a regra do 
trapézio. Tem-se agora 


[=> (Ea) + (a) + E + (Fay + 
+Eb)+E” (3) 


ou desprezando E' + E”, 


I=Lta)+ta)+ Eb) (4) 
2 2 


com um erro bastante menor do que o primeiro. 
Se dividirmos o intervalo (a,b) em muitas par- 
tes iguais — n partes iguais, de comprimento h — 
obtém-se um valor cuja aproximação cresce 
com n; para funções suficientemente regulares, 
o erro pode tornar-se tão pequeno quanto se 
quiser; quando o sentido da concavidade da 
curva mude, há uma compensação de áreas, em 
excesso e em defeito que apressa a redução do 
erro. Chega-se assim à regra do trapézio para o 
cálculo dos integrais definidos. 


p= (Ly ton) (5) 


onde os yy são os valores da função nos pontos 
de divisão do intervalo. A fórmula mostra que, 
conhecidos esses valores (a função pode ser, por 
exemplo, determinada experimentalmente), as 
quatro operações bastam para o cálculo do inte- 
gral. Quando a função seja dada numa tabela, o 
cálculo dos valores não tabelados pode fazer-se 
por fórmulas de interpolação aritméticas que só 
envolvem, também, as quatro operações. 


REGISTOS, MEMÓRIAS. NÚMERO E COM- 
PRIMENTO. O computador ABCD pode consi- 
derar-se à primeira aproximação como um con- 
junto de registos ou locais, dos quais os mais 
importantes se chamam memórias. Um registo é 
(Figura 2) um lugar onde se formam fieiras de 
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Fig. 2 — Um registo ou memória 


dígitos, sequências de realizações materiais de 
dígitos, e se pode imaginar (vejam-se as janelas 
das máquinas de calcular manuais) como fieiras 
de algarismos árabes. O acumulador duma má- 
quina de secretária é um registo. Além de dígi- 
tos também se realizam sinais, como se vê na 
Figura 2. O conteúdo de uma memória chama-se 
palavra ou termo. As memórias estão numeradas e 
o número de cada memória é o seu endereço. 
Representemos por 123) o conteúdo da memó- 
ria 123, por x o da memória X; para significar 
o endereço da memória que contém o número 123 
escreveremos (123) ou loc 123 (de locus, lugar). 

Para o utilizador, as primeiras características 
de um computador, relativas à sua capacidade, são 
o número de memórias e o comprimento das 
palavras. O computador ABCD tem 1000 memó- 
rias, cujos endereços são 000, 001,... 999; a 
palavra é de 13 dígitos decimais, e sinal. Os re- 
gistos de memória, ou memórias, da máquina 
formam a sua memória. 

Quando se queira representar o conteúdo de 
partes da memória do computador escrever-se-á, 
por exemplo 


567) + 1124125125568 
568) + 1123124124569 


571) + 1124126127678 


ou, por motivos que imediatamente se verão, 


567) + 1 124 125 125 568 
568) + 1 123 124 124 569 


571) + 1 124 126 127 678 


Além das memórias onde se inscrevem os 
dados, o programa de cálculo e os resultados, o 


computador deverá ter circuitos de interpretação 
do programa, de execução de cálculos e outras 
operações, de contrôle geral, que estão, por 
enquanto, fora dos nossos interesses; os conteú- 
dos da memória em épocas sucessivas (segundo 
uma cadência característica do computador) cons- 
tituem os instantâneos dos passos sucessivos no 
processo de cálculo. Depois do número de me- 
mórias que um computador tem e de quantos 
dígitos tem cada palavra, importa saber o sis- 
tema de numeração que a máquina usa. Na me- 
mória do computador ABCD há em cada época 
1000 palavras. Cada palavra é constituída por 13 
dígitos decimais (decits); além disso há uma 
14.2 posição à esquerda onde figura um sinal 
algébrico. A aritmética que utiliza é a decimal. 

Uma palavra pode ser um número ou o código 
de uma instrução. Como se entenderá melhor mais 
adiante, esta ambiguidade sobre que se fundam 
muitas das virtualidades dos computadores é des- 
feita pela máquina, que trabalha a dois tempos. 
Num dos tempos, o tempo de instrução, a máquina 
vai a uma memória procurar que fazer e toma o 
seu conteúdo como um código que decifra nos 
circuitos adequados; no outro tempo, o tempo de 
execução, a máquina, procurando numa memória 
um número, interpreta como tal o que lá encontra, 

Assim 

567) + 1124125125568 


quer dizer que a memória 567 contém o número 
0,1124125125568 


ou a instrução «some o conteúdo da memória 124 com 
o conteúdo da memória 125 e ponha o total na memó- 
ria 125; vá buscar à memória 568 a instrução que 
cumprirá a seguir». Neste último caso preferirá 
escrever-se 


567) 1 124 125 125 568 


Como se disse, a ambiguidade não subsiste 
no processo de cálculo; se a máquina, tendo aca- 
bado a execução de uma ordem, vai ao registo 
567 buscar uma instrução, analisa o conteúdo 
do registo supondo-o uma instrução ; se, na exe- 
cução de uma ordem, a máquina vai buscar ao 
registo 567 um número, interpreta o termo como 
sendo o número que procura. No seguimento 
omitiremos, como aqui, o sinal nas instruções. 

As memórias retêm as palavras que nelas se 
formam emitindo-as para os circuitos aritmé- 
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ticos ou outros registos no decurso dos cálculos 
quando a máquina interpreta e executa uma ins- 
trução que comanda essa emissão. A substitui- 
ção do conteúdo de uma memória é realizada 
por instrução, também. 

Quando se comanda a inscrição dum termo 
numa memória o seu conteúdo é previamente 
rasurado; mas quando se comanda uma opera- 
ção em que intervém o conteúdo duma memória, 
este conteúdo permanece. A maneira, portanto, 
de «limpar» uma memória é comandar a inscri- 
ção do zero. 

Dos registos que não são memórias o mais 
importante é o acumulador, que se representa ge- 
ralmente por A como o seu conteúdo por a. 
O acumulador é um registo aritmético onde se 
formam resultados de operações. No computa- 
dor ABCD o acumulador não é endereçável e não 
tem que interessar-nos como utilizadores. 


ESTRUTURA DAS INSTRUÇÕES. Posto isto, 
que informação deve o computador encontrar 
codificada numa instrução para executar uma 
operação aritmética? Suponha-se que o compu- 
tador está construído para executar as quatro ope- 
rações, interessando dois dados, (a e b) para obter 
um resultado (c). Cada instrução deve então 
indicar a operação, os endereços dos dois dados 
(A e B) e o endereço da memória onde o resul- 
tado será inscrito (C). 

Veja-se que, uma vez cumprida a instrução 
corrente, o computador tem de saber onde en- 
contrar a que deverá cumprir seguidamente. No 
computador ABCD o endereço da memória onde 
o cumputador vai procurar o que fazer, depois 
de cumprir a instrução, é indicado nesta instru- 
ção por D. 

Quando uma palavra é uma instrução os seus 
13 dígitos são utilizados na codificação da infor- 
mação agrupados da esquerda para a direita 
como indica a decomposição 


3=1+3+3+3+3 


O primeiro decit indica a operação segundo 
a correspondência 


1 adição 

2 subtracção 

3 multiplicação 
4 divisão 
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Os três seguintes compõem o endereço A do 
dado a; os três seguintes, o endereço B do 
segundo dado; os outros três, o endereço € do 
resultado. Finalmente, os últimos três compõem 
o endereço da memória onde se encontra a ins- 
trução seguinte. Este quadro resume o que acaba 
de dizer-se: 


CÓDIGO ACÇÃO 

1 ABCD a+E —s € = D 
2 ABCD ED me E; cel do D 
3 ABCD aseb -—s. E; ef de D 
4 ABCD at — E =| ke D 


EXEMPLO. Assim, suponhamos que em certa 
época o contrôle está no local 567, isto é, que a 
máquina vai explorar 567 em busca de uma 
instrução a cumprir. Imaginemos que a memó- 
ria apresenta o seguinte aspecto 


123) 2 
124) 1 
125) 0 
567) 1 124 125 125 568 
568) 1 123 124 124 569 
569) 1 124 125 126 570 
570) 1 123 124 124 571 

1 124 126 127 678 


571) 


O conteúdo de 567, analisado pela máquina 
como uma instrução, é interpretado como a 
ordem de adicionar o conteúdo de 124 ao de 
125, de colocar o resultado em 125 e de pro- 


curar em 568 a instrução seguinte: 


124) + 125) — 125, LU. 568 


A máquina vai pois somar 1 a O e pôr o 
resultado, que é 1, em 125. 

Os registos de memória, lidos quando o seu con- 
teúdo intervém numa operação, conservam-no, 
como se disse, até que recebam a ordem de ins- 
crição do conteúdo de outra memória, que pode 
ser de rasar (correspondente a uma inscrição de 
zero), ou de uma leitura do exterior. Portanto, 
executada a instrução 567), a única alteração no 
panorama da memória é agora que 125) é 1, em 
vez de O. Tem-se, pois, 


123) 2 124) 1 125) 1 


Passando à instrução 568), o computador in- 
terpreta-a como 


123) + 124) — 124, LN 569 
e o panorama é, no fim da sua execução, 


123) 2 124) 3 125 1 


Continuando, executada 569), tem-se 


123) 2 124) 3 125) É 
126) 4 
no fim de 570) 
123) 2 124) 5 125) 1 
126) 4 
e no fim de 571) 
123) 2 124) 5 125) 1 
126) 4 
127) 9 


Observe-se que as instruções 568) e 570) 
colocaram em 124 os números 3 e 5; uma ins- 
trução 572) idêntica faria aparecer 7. O pro- 
grama pode pois ser prolongado por forma a 
gerar em 124 a sucessão dos números ímpares, 
enquanto houver espaço para armazenar as ins- 
truções correspondentes. Por outro lado, as ins- 
truções 567), 569) e 571) fizeram aparecer a 
partir do local 125 os quadrados 1, 4, 9. Como se 
sabe, desde Pitágoras, as somas parciais obtidas 
na série 


1+-3--p5+7+H.... 


constituem a sucessão dos quadrados perfeitos, 
como recorda a Figura 3 


++ ++ 
+ tt + 6 É d+ 

+ + oo +oo io + 
o o + oo + ooo «+ 


Fig. 3 — Cálculo de quadrados 


Exercício 1. Prolongar o programa anterior por 
forma a calcular os quadrados dos primeiros 10 
inteiros. 


PROGRAMAÇÃO DE OPÇÕES. A primeira 
instrução, que ultrapassa a automatização linear 
de uma sucessão de operações aritméticas, abrindo 


assim promissoras possibilidades é a da progra- 
mação de opções: A instrução é (!) 


S ABC =É lo C sss 
ac: Leno) sss 


a>b 
a< b 


Trata-se de uma instrução de contrôle: a máquina 
recebe a ordem de ir procurar, analisar e cum- 
prir uma outra instrução. Esta outra instrução pode 
ser uma de duas, ou a instrução contida em € 
ou a instrução contida em D, e o que decide 
qual das instruções vai ser cumprida é a com- 
paração dos conteúdos dos registos A e B. Se o 
número em A for maior que o número em B, a 
máquina transfere o contrôle para C; se não, o 
contrôle é transferido para D. Ora em Ae B 
podem estar resultados de cálculos anteriores o 
que dá à máquina a aparência de uma autono- 
mia, que de facto não é mais do que a execução 
do que se contém na instrução 5. 

Esta instrução é uma instrução de transferên- 
cia (ou salto) condicional. Quando houver que 
transferir o contrôle para a memória Z, incon- 
dicionalmente, bastará usar a instrução 


52222 


Exercicio 2. Verificar que a instrução 


52222 
é uma instrução de salto incondicional para Z. 


Outra instrução de contrôle é a ordem de 
paragem 
OABCD 


onde os endereços são mudos e se pode repre- 
sentar sempre por uma sucessão de 13 zeros. 
Suponha-se o pormenor do programa 


789) 5 100 200 790 791 
790) O 


Quando na execução o controle passasse para 
789, e o conteúdo de 200 não atingisse o número 
guardado em 100, o cálculo pararia. No caso 
contrário, o contrôle passaria para 791 onde 
poderá principiar um novo cálculo. Este cálculo 
só se efectuará, portanto, se o número em 200 
(que pode ser resultado dum cálculo anterior) 


(!) A abreviatura sss quer dizer «se e só se» ou «sem- 
pre que e só quando». 
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não for menor que o número em 100 (que pode 
ser uma constante introduzida inicialmente. 


LEITURA E ESCRITA. As ordens de input- 
output, de cópia (para dentro e para fora do 
computador), ou de leitura ou de escrita, são 


6ABCD — AB MN. D 
7ABCD | ABL LD 


A instrução 6 manda ler para as memórias de 
A a B, incluídas ambas, e transfere o contrôle 
para D. A instrução 7 manda escrever das memó- 
rias de A a B, incluídas ambas, e transfere o 
contrôle para D. 

Vamos supor que o computador lê fita perfu- 
rada e escreve perfurando fita. A fita pode ser 
de quatro fileiras, representando-se em binário 
os dígitos, como na Figura 4, e os sinais por 
duas composições disponíveis. 


0 5 + 2. + 
1 + 6 + +. 
2 + 7 a po + + 
3 + + 8 + 
4 - . o ap +- 
POSITIVO +44 ++ 
NEGATIVO +++. 


Fig. 4 — Código de perfuração 


O computador, quando decifra uma ordem de 
leitura, activa o leitor de fita e pára se não hou- 
ver fita colocada; havendo fita, esta corre até ao 
primeiro sinal de + ou —. Este sinal e os dígi- 
tos seguintes (até ao próximo sinal e no máximo 
de 13) são colocados na primeira memória (A) 
indicada na instrução, o próximo sinal e os seus 
dígitos ma memória seguinte, e sucessivamente 
até se ter preenchido a segunda memória (B), 
indicada na instrução; o contrôle passa então 
para D. O computador pára de ler se a fita se 
esgota, o que permite ler um número ignorado 
de palavras. 

Note-se que, além da leitura comandada por 
instruções de programa, tem que haver leitura 
comandada exteriormente — pelo menos para fazer 
entrar um programa de leitura. O computador tem, 
num painel de comando, um botão de leitura L 
e uma série de botões onde se inscreve, com os 
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dígitos decimais codificados em binário (diz-se, 
em DCB), o endereço E da primeira memória a 
preencher, como se mostra na Figura 5 (onde 
se tem representado o endereço 567, premidos 
os botões t+). Além do botão L, o painel tem 
ainda o botão I que serve para principiar a exe- 
cução do programa cuja primeira instrução esteja 
na memória E; começado o cálculo, outros sub- 
programas podem ser cumpridos a partir de 
ordens de transferência para os locais onde a 
sua inscrição se inicia. 


. + k + 
++. 
++ + 
L I 

Leitura Início 


Fig. 5 — Pormenor do painel de contrôle do computador 


OPERAÇÃO DO SISTEMA. Por muito impor- 
tante que seja a parte que o computador tem na 
actividade de um centro de cálculo, outras enti- 
dades são igualmente indipensáveis numa ou 
noutra estação do processo de cálculo. O problema 
que o computador resolve começou por ser, por 
exemplo, um problema de física. A análise mate- 
mática, mais ou menos complexa, leva à sua 
formulação como um problema de cálculo ; nesta 
altura, quere-se determinar um número ou um 
sistema de números. Poderá dizer-se que se trata 
sempre de resolver uma equação ou um sistema 
de equações, mas esta generalidade é inútil; 
convém antes ver se o problema cabe num dos 
seguintes grupos, por exemplo, ou se é suscepti- 
vel de ser decomposto em problemas que caibam, 
o que acontecerá quase sempre: 


1. Cálculos de valores de funções dadas por 
fórmulas ou tábuas (empíricas ou calculadas). 

2. Solução de sistemas de equações lineares; 
determinação de valores próprios e vectores pró- 
prios de matrizes; outros cálculos ligados à álge- 
bra matricial. 

3. Solução de equações não lineares, ou de 
sistemas destas equações. 

4, Cálculo de valores de funções dadas por 
equações diferenciais (ordinárias ou parciais) e 
condições nos limites (fronteiras no espaço ou 
épocas no tempo). 


5. Cálculos de derivadas de funções em pontos, 
ou de integrais em intervalos, dadas as funções 
por tábuas. 


6. Cálculos de valores de funções dadas por 
equações integrais. 

7. Cálculos envolvendo compensação e filtra- 
gem, como os de mínimos quadrados e de análise 
harmónica, clássica ou generalizada. 

8. Cálculos de extremos, como os da progra- 
mação linear. 


Analisado o problema de cálculo, há que es- 
colher os algoritmos, as técnicas de aritmetização 
aplicáveis em cada ocasião e ainda que estabelecer 
a concatenação das partes, o correr do processo. 
Aquela escolha põe problemas de análise numé- 
rica; esta concatenação constitui a ordenação e 
faz-se, as mais das vezes, com organigramas. A 
aritmetização e a ordenação do problema de 
cálculo é a programação própriamente dita. Segue-se 
depois a codificação das rotinas parciais e a sua 
integração no programa, implicando a codificação 
das articulações. A aritmetização, a ordenação, 
a codificação e a integração constituem a progra- 
mação no sentido lato de que sai o programa. Os 
computadores modernos são por vezes descritos 
como máquinas aritméticas de programa interno ; 
quer-se com isto sublinhar que a máquina obe- 
dece durante a execução do cálculo a um programa 
que está inscrito nos seus registos. A perfuração 
do programa e a sua leitura é portanto o passo 
que se sucede à integração no programa (com 
papel e lápis) das rotinas codificadas. 

Quanto aos dados, nada impede que estejam 
incluídos no programa. Repare-se que num pro- 
grama de cálculo há sempre constantes a consi- 
derar. Quando se queira, por exemplo, calcular 
valores de um polinómio, o seu e os coeficientes, 
bem como os valores do argumento, são cons- 
tantes a dar ao computador. Agora, várias situa- 
ções podem ocorrer: ou se quer um programa 
para uma ocasião ; ou se prevêem várias ocasiões 
em que o polinómio é o mesmo e diversos os 
valores do argumento; ou se quer um programa 
para qualquer polinómio e quaisquer valores 
(a serem dados oportunamente) ; etc. 

Como, todavia, o trabalho de programação é 
um esforço considerável, um encargo pesado e 
até uma tarefa difícil, todos os centros tendem 
a constituir uma biblioteca de programas capaz de 
resolver os problemas que se repetem. Em cada 


caso o programa está perfurado e pronto a ser 
lido; só resta perfurar os dados. Na verdade, os 
utilizadores de cada máquina tendem a organi- 
zar-se para pôrem em comum os seus programas, 
e uma máquina vale não só pelas suas caracte- 
rísticas como pelas instituições que as utilizam. 
Por exemplo, o Pegasus uns anos atrás e o 
NATELL803 agora, só porque foram escolhidos 
por várias universidades inglesas, onde existe 
tradição valiosa de trabalho sério, trazem con- 
sigo bibliotecas de programas muito valiosas 
para países pobres. Há ainda um segundo mótivo 
para separar os programas dos dados, motivo 
que se refere ao espiolhamento dos programas, se- 
guinte à perfuração. Não se tratando de uma 
actividade directa de cálculo, pois que se trata 
da pesquisa dos erros de programação, é uma 
actividade considerável na realidade. Para auxi- 
liar o espiolhamento dum programa preparam-se 
os dados de um exemplo, cujos diversos passos 
se possam obter convenientemente sem o uso do 
computador; a execução completa do cálculo 
simplificado que serve de exemplo precede, natu- 
ralmente, a operação do programa com os dados 
reais. Os dados reais são neste caso perfurados 
independentemente, embora eles e o programa 
possam ser fâcilmente reproduzidos numa só fita 
quando se queira. 

Tudo isto leva à separação de duas ordens de 
trabalhos junto de um computador: a programa- 
ção para o computador e a operação do computador. 

a primeira ordem de trabalhos que temos 
vindo a descrever; nela a leitura do programa é 
seguida da leitura dos dados do exemplo e da 
tentativa de cálculo. Entra-se assim no ensaio do 
programa, que termina logo se a tentativa é bem 
sucedida, o que só raramente acontece. Quando 
não é assim, o programador procura fazer o 
diagnóstico das causas do erro, o que o leva a 
um momento anterior do processo de programa- 
ção para introduzir correcções cujas conseguên- 
cias deverá seguir dali para a frente. Em geral 
só foi descoberto um erro e o programa tem vá- 
rios, sucedendo-se os retornos até que o exem- 
plo é levado do princípio ao fim correctamente. 


CÁLCULO DE INSTRUÇÕES. Como se disse, 
o número de registos de memória de um compu- 
tador (a capacidade da sua memória) é uma ca- 
racterística que, tanto como as velocidades a que 
se cumprem as instruções elementares, condi- 
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ciona a velocidade efectiva do cálculo. Na ver- 
dade, quando se diz que um computador faz 1000 
adições por segundo deveria talvez dizer-se que 
adiciona à cadência de 1000 operações por se- 
gundo. Estando as parcelas memorizadas no iní- 
cio do cálculo, um computador com uns milhares 
de memórias só pode executar uns milhares de 
adições de uma assentada. Se, pelo contrário, as 
parcelas foram introduzidas oportunamente, é a 
duração da leitura que reduz drâsticamente aquela 
celeridade. 

Num sentido largo, e todavia justo, pode 
mesmo considerar-se que o tempo de programa- 
ção é ainda uma parte do tempo de cálculo ; ora, 
sendo um cérebro electrónico perfeitamente estúpido, 
isto é, completamente incapaz de fazer o que não 
esteja programado, a programação é uma tarefa 
delicada, difícil e morosa. Acrescendo que dos 
computadores se exige a resolução de problemas 
que com outros meios estão fora do alcance, um 
ou dois anos de vários programadores vão às 
vezes no mesmo programa. Mesmo no tempo de 
execução dos cálculos, o tempo consumido no cum- 
primento das instruções que não são de cálculo, 
mas de ordenação do processo, é com frequência 
uma parte importante. 

Voltemos, porém, ao caso simples e raro em 
que só a capacidade e a velocidade da máquina 
contam (como no caso da tábua de quadrados). 
Como vimos, cada quadrado é obtido por duas 
instruções, uma que calcula o novo ímpar, outra 
que soma este ímpar ao quadrado anterior. Con- 
tando com o local onde se guarda o quadrado, 
tem-se um total de 3 registos ocupados no seu 
cálculo; o programa ocupa o dobro da memória 
onde se guardam os resultados. Assim, no com- 
putador de 1000 memórias, que estamos imagi- 
nando, não poderão calcular-se de uma só vez 
muito mais de 300 quadrados com aquele pro- 
grama. 

Já se insistiu na perfeita ambiguidade das pa- 
lavras que faz com que o conteúdo de um local 
seja uma instrução ou um número conforme o 
tempo em que a máquina o explora. Vamos agora 
ver como essa ambiguidade, permitindo o cálculo 
de instruções possibilita consideráveis encurtamen- 
tos nos programas onde as mesmas instruções 
se repetem com simples modificações de ende- 
Teços e OS programas constituem, a menos destas 
modificações, a repetição de um mesmo ciclo de 
cálculos. 
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Exercício 3. Estudar o processo de cálculo que 
começa com o contrôle no registo 111 sendo 


My LH. 7a Mi 


Seja outra vez o programa 


567) 1 124 125 125 508 
568) 1 123 124 124 569 
569) 1 124 125 126 570 
570) 1 123 124 124 5711 


A instrução 569) difere essencialmente da ins- 
trução 567) em que o terceiro endereço é 126 
em vez de 125; seguidamente, a instrução 570) 
não difere essencialmente de 568). Se, portanto, 
em vez de se passar ao cumprimento de 569), 
se puder substituir 125 por 126 em 567) e cum- 
prir de novo esta instrução, tem-se o mesmo 
resultado. Ora, visto que o conteúdo de 567 se 
pode considerar um número, se for cumprida a 
instrução 


E) 1 567 (1000) 567 507 


passa a ter-se, efectivamente 


56T) 1 124 125 126 568 
568) 1 123 124 124 569 


Pondo a instrução E) em 569 será inicialmente 


567) 1 124 195 12 568 
568) 1 123 124 124 569 
569) 1 567 (1000) 567 567 


que prolongamos representando por 567)' o con- 
teúdo de 567 depois da primeira modificação de 
endereço, etc. 


567) 1 124 125 126 568 
568) 1 123 124 124 569 
569) 1 567 (1000) 567 567 
567)” 1 124 195 127 568 


Como se reconhece, embora 567)' seja exacta- 
mente o que convém, o mesmo não se passa 
com 567)” que deveria ser 


567)” 1 124 126 127 568 


O que acaba de acontecer é característico do 
que se passa na programação dos pormenores 


de cálculos (pelo menos para programadores de 
temperamento mais vulgar): levantando um pro- 
blema, aparece uma ideia que quase o resolve, 
isto é, que o resolve desde que se façam ajusta- 
mentos; o programa fica mais complicado do 
que fazia prever a sugestão inicial e esta chega 
a ficar irreconhecível na sucessão final das ins- 
truções. Assim, a ordenação dos cálculos das 
partes faz-se de fora para dentro, da considera- 
ção do problema geral para a definição dessas 
partes; a codificação, pelo contrário, parte de 
pontos bem definidos na região do cálculo que 
se considera e segue de dentro para fora, intro- 
duzindo a pouco as complicações que vão resol- 
vendo problemas residuais. 
É claro que se for 


569) 1 567 (1001000) 567 567 


o programa é o que convém, excepto da primeira 
vez em que a instrução é cumprida. A solução 
da dificuldade estará em alterar as instruções 
iniciais, fazer o que se chama a preparação do 
ciclo. Basta para isso, como fácilmente se reco- 
nhece, que se parta do estado inicial. 


123) 2 124) 3 125) 1 

e se mande 
567) J 124 125 126 568 
568) 1 123 124 124 569 
569) 1 567 (1001000) 567 567 
Teremos 
567) 1 124 126 127 568 


que encontra 


123) 2 124) 5 125) 1 
126) 4 

e deixa 
123) 2 124) 5 125) 1 
126) 4 
127) 9 


O programa vai executar o cálculo de quadra- 
dos enquanto este cálculo não violar algum por- 
menor de construção, ou o programa se não 
alterar a si próprio pela repetição não limitada 
da modificação do endereço. É preciso programar 
um fim de cálculo. 


Exercício 4 Descrever o desenrolar do programa 
anterior. 


CONTAGEM E ENSAIO. Quando um ciclo 
deva repetir-se n vezes terá que programar-se: 
i) a contagem do número de repetições concluí- 
das; ii) o ensaio que no fim de cada iteração de- 
terminará se n foi atingido ou não; iii) a trans- 
ferência (para o princípio do ciclo enquanto se 
não tenha atingido n, e para o seguimento 
quando se tenha atingido n). 

O contador é uma memória K, que no início 
do cálculo contém 2 (por motivos que se verá) 
e a cujo conteúdo se soma uma unidade por 
cada repetição. A instrução 


O Lt RE aE 

soma ao conteúdo de K uma unidade sempre 
que seja preciso, transferindo depois o contrôle 
para X. 

O ensaio compara k com n transferindo o con- 
trôle para um princípio de ciclo quando n>k 
e para o seguimento (que, suponhamos, começa 
em 5) quando n==k, Recorde-se que a instrução 


5SABCD 
comanda a acção 
a Vs C 5 25+5 
cade las sss a>b 


A instrução do teste é, então, 


TT SNK Ss 

A transferência para KI significa que se vai 
aumentar o número de repetições da unidade 
correspondente à repetição que se inicia imedia- 
tamente. O último endereço em KI deve pois 
ser 567. 

Pode agora fixar-se os endereços simbólicos pondo, 
por exemplo, 


T = 568 
KI = 569 


N= 121 
K = 122 


o que arrasta ajustamentos nos restantes ende- 
reços. O programa para o cálculo e perfuração 
de n quadrados é, agora, 


567) 1 124 125 126 568 
568) 5 121 122 569 572 
569) 1 (1) 122 122 570 
570) 1 123 124 124 571 
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S21) 1 567 (1001000) 567 567 
572) 7 Z%25 125 +-n 573 
573) O 


Os computadores ABC 


ECONOMIA E SIMPLICIDADE. Quando se 
revéem os planos de um computador com o 
cbjectivo de se introduzirem melhoramentos, são 
vários os objectivos que se perseguem, vários os 
critérios por que se aferem esses melhoramentos 
e diversos os motivos que os tornam possíveis. 
Assim, os progressos da tecnologia podem tor- 
nar possível a realização de circuitos menos dis- 
pendiosos, mais seguros, mais sóbrios no que 
consomem e mais resistentes ao uso. Também, 
por outro lado, certos progressos de concepção 
podem trazer economias nos custos ou simplifi- 
cações na programação, dois objectivos sempre 
muito importante. 

Deste ponto de vista da programação (pois 
que ainda nos não ocupámos de como são feitos 
os computadores) façamos umas observações ao 
anterior. Em primeiro lugar, é patente que a ins- 
trução seguinte está quase sempre no registo se- 
guinte ao da instrução presente, pelo que ocorre 
arranjar as coisas por forma a só escrever o seu 
endereço, durante a programação, quando não 
for assim. O computador é uma máquina onde 
se realizam convenções; pode estar de tal modo 
construído que, cumprida uma instrução, procure 
a nova instrução no registo seguinte — salvo 
quando seja explicitamente instruído a fazer 
outra coisa. Nestas condições cada instrução 
poderá dispensar os últimos três dígitos do re- 
gisto, utilizando só 10. E, como 10 dígitos dão 
uma precisão suficiente, faz-se uma grande eco- 
nomia reduzindo o tamanho da palavra a 10. 
Assim na sua segunda versão o computador 
ABCD, será um computador ABC. 

A única instrução que precisa ser alterada é a 
transferência condicional 


SABCD 


que passará a ser 


SABC 


com a seguinte acção: quando seja a»b o con- 
trôle passa a C; quando a<b o contrôle con- 
tinua em sequência. 

Diz-se de um computador assim, que tem o 
contrôle em seguência. 
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Todas as vezes que se tiver que dar um salto 
incondicional usa-se mais uma instrução com 
perda de espaço e tempo, mas na grande maio- 
ria dos passos de cálculo poupa-se material. 
Quanto ao tempo, depende do tipo de memória 
do computador porque nem todas as memórias 
são adequadas ao contrôle em sequência. 

Há uma meia dúzia de anos, o tambor magné- 
tico constituía a memória principal de muitos 
computadores. Assim, por exemplo, a memória 
principal da ZEBRA, do equipamento LNEC de 
1959 a 1963, é constituída por um tambor mag- 
nético cilíndrico, rodando à considerável veloci- 
dade de 100 rps por forma a fazer passar rapi- 
damente as suas 256 pistas em frente de outras 
tantas cabeças de leitura. Cada pista tem assim 
em série 32 memórias numeradas seguidamente. 

Se o computador ABC tiver um tambor com 
estas características, cada volta dura 10 milise- 
gundos e só 10 milisegundos depois de lida com- 
pletamente uma instrução pode começar a ser 
lida a instrução de seguência (pensemos que a 
interpretação e execução consomem tempo). Os 
tempos de interpretação e execução de uma or- 
dem sendo assás variáveis na escala dos milise- 
gundos, p>derá contar-se, em média, com 5 mi- 
lisegundos de espera pelo registo seguinte na 
cabeça de leitura, o que não é insignificante em 
face dos 312 microsegundos que leva a ser lida 
uma palavra. 

Como, por outro lado, a execução envolve em 
geral dados inscritos noutras memórias, entende- 
-se que numa máquina de tambor magnético a 
programação em sequência mantenha inactivas 
as cabeças de leitura numa parte considerável 
do tempo. Esta perda de tempo é agravada, rela- 
tivamente, pela circunstância de ser motivada 
pelas características do movimento do tambor; as 
velocidades muito maiores dos processos electró- 
nicos acomodam-se mal com a existência de ele- 
mentos móveis, dos quais exigem performances 
extraordinárias que custam caro. 

As máquinas de tambor magnético têm uma 
programação chamada óptima que não é sequen- 
cial. Nela se leva em linha de conta os tempos 
decorridos antes de se fixar o endereço da me- 
mória onde deve estar, oportunamente, o dado 
ou a instrução de que então se precise. 

O contrôle em sequência que estamos a ima- 
ginar para o computador ABC seria pouco efi- 
ciente numa máquina de tambor, pelo que vamos 


supor que ele tem uma memória de núcleos 
magnéticos em que todos os registos, que se 
atingem no mesmo tempo, são igualmente aces- 
síveis. Esse tempo é no National-Elliott 803, por 
exemplo, de 288 microsegundos, o que corres- 
ponde a economia apreciável no tempo de espera, 
ou a programações mais simples, como se queira. 

Imaginamos assim transformado o computa- 
dor ABCD, que era uma máquina com pala- 
vras de 13 dígitos e instruções de 3 + 1 ende- 
reços, no computador ABC, que é uma máquina 
com palavras de 10 dígitos e instruções de 3 ende- 
reços. Ambas as máquinas têm 1000 memórias. 
A primeira, num tambor magnético; a segunda, 
em núcleos magnéticos. Vamos também supor 
modificado o computador por forma a aumentar 
o número de instruções ao seu código, facilitando 
a programação e a economia de espaço na me- 
mória. Teremos, naturalmente, de aumentar e 
complicar os circuitos de análise e execução, e 
portanto, de encarecer o computador ; este agra- 
vamento de preço será o custo de uma maior 
eficácia, num sentido muito geral em que se con- 
sidere a facilidade de programação e o tamanho 
dos cálculos. Haverá ainda que ter em conta os 
tempos, que podem diminuir pelos progressos 
tecnológicos mas aumentar com a simplificação 
do código; sem esquecer que os tempos de pro- 
gramação também pesam. 

A introdução das novas instruções vai ser es- 
tudada a propósito da programação do cálculo 
de polinómios. 


APROXIMAÇÃO POR POLINÔMIOS. O cál- 
culo de funções de uma variável utiliza geral- 
mente, de uma maneira ou outra, a substituição 
por um polinómio que difere da função, para os 
mesmos argumentos, de quantidades que se des- 
prezam. Quando sejam dados os valores da fun- 
ção e suas derivadas na origem, a fórmula de 
Maclaurin 


FO), 


1! 


f (x) = f(O) + fo 


pe tg, x" + Rn 


- EM (O) 
+ ——— 
2! n! 


permite pôr 


f(x)=aotax+ax'+H...+tanx" 


quando, para os valores que interessam seja des- 
prezível o resto Rn. É então 
EMO) 
TO 
Outros polinómios que aproximam convenien- 
temente funções são os polinómios interpoladores 
que se fazem coincidir com ela em vários pon- 
tos, quando a função seja assim dada. Podem 
obter-se a partir das tabelas de diferenças da 
função, por várias fórmulas ou pela combinação 
linear de polinómios de várias famílias; a fór- 
mula de Maclaurin descreve-se desta posição 
como uma combinação dos polinómios de base 


1xx]...x" 


aproveitando a informação sobre o comporta- 
mento profundo da função na origem fornecido 
pelas derivadas. 

Da mesma maneira, a fórmula de Lagrange 
pode dizer-se que dá o polinómio que toma va- 
lores dados em certos pontos como combinação 
linear dos polinómios que se anulam em todos 
excepto num (onde tomam o valor 1). 

Uma família de polinómios convenientes para 
este fim, quando só interessem valores do inter- 
valo (0,1), é a dos polinómios de Tchebychef 
(Tn). Na Figura 6 tem-se os coeficientes de al 
guns polinómios, supostos ordenados por potên- 
cias decrescentes, em face dos respectivos graus n. 


n 
Õ 1 

1 2 — | 

2 o — 8 1 

5) 92  —48 18 e | 

4 128 —256 160 —32 l 

5 62 — 1280 1120 —400 50 -—1 

O 2048 — 6144 6912 —3584 840 —72 1 


Fig. 6 — Coeficientes dos polinómios T, 


Os polinómios de Tchebychef podem obter-se 
das tórmulas 
T. (x) = cosn9 cosg = 2% —1 
onde se vê que tem de ser 0< x< 1. 


Exercício 5. Verificar, partindo destas relações, 
quatro linhas do quadro anterior. 
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CÁLCULO DE POLINÓMIOS. A programa- 
ção do cálculo de polinómios pode fazer-se com 
variados objectivos, como se poderá fazer ideia 
das situações que se distinguem a seguir. 


[ — Polinómio do 4.º grau dado. Valores dados 
do argumento devendo ser lidos quando 
precisos. Impressão imediata dos resulta- 
tados. 

Il —O mesmo, com leitura inicial de todos os 
valores do argumento e escrita, no fim de 
todos os valores calculados. 

[HI— O mesmo, com os valores do argumento 
em progressão. Por exemplo, 


x== 0,00; (0,01); 1,00 


IV — Polinómio, cujo grau (não superior a 10) 
será dado com os coeficientes. Valores do 
argumento (em progressão aritmética) a 
serem dados). Leitura e escrita na ocasião. 

V —Polinómio do 4.º grau dado. Valores do 
parâmetro, calculados num programa prin- 
cipal, em número dado. 

VI — Polinómio de grau não superior a 10 a ser 
calculado num programa principal. Valor 
do parâmetro também calculado no pro- 
grama. 

VII—100 polinómios de grau 10, dados. 

O processo de cálculo utilizado é o processo 
das cascas que se funda na identidade 


ax"taxtit...+tanmxta = 
= (((ax+ta)x+ta)x+...)x+Han 


onde se reconhece um ciclo de operações, como 
convém. Seja P o endereço da memória onde se 
quer formar o polinómio; A,,..., An, as memó- 
rias onde estão os coeficientes, X o endereço 
de x. Suponhamos que se coloca a, em P no 
início do cálculo. Então as instruções 


M) MUL P X P 
Si) SOM P A P 


constituem um ciclo que repetido n vezes, para 
i= 1, 2,...n, nos coloca em P o valor do poli- 
nómio. 

Observe-se que utilizou uma simbologia (cuja 
clareza dispensa palavras) para representar a 
adição e a multiplificação. Este código mnemónico 
do código da máquina tem sobretudo vantagens 
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quando se aumenta o número de instruções e o 
código da máquina deixa de se reter fácilmente. 


PROGRAMAÇÃO BANAL DO CÁLCULO 
DE POLINÓMIOS. Seja o cálculo do polinómio 


P(x)=128x*— 256xº + 160x?— 32x +1 


para x = 0,1. 


Leituras. Os coeficientes ao ,-.., a; serão lidos 
para as memórias 301 a 305 por uma ordem 


L1) LES Sol 305 
o de x para 201, por 
L2) LEI 201 201 


e o programa para 101, e seguintes, por uma 
ordem 


L3) LEI 101 FIMM 

onde se representa por FIMM o endereço da 
ordem de paragem que será cumprida quando o 
cálculo estiver terminado. (Adoptamos a prática 
de utilizar endereços simbólicos, que nunca serão 
de três letras para não confundir com os códigos 
simbólicos de operações, que serão sempre de 
três letras). 


Programa. Vamos começar pelo programa arit- 
mético a partir dos dados 


x ão a, a, a; a, 
501) 
101) MUL 301 201 501 ax 


102) SOM 501 302 501 aox+x, 
103) MUL 501 201 501 (ax+a)x 
104) SOM 501 303 501 


(ax a,)x + a; 
105) MUL 501 201 501 ... 
106) SOM 501 304 501 
107) MUL 501 201 501 ... 
108) SOM 501 305 501 P(x) 


Escrita e paragem. Fazem-se pelas instruções 


109) IMP 501 501 000 


110) PAR 


MAGUE 


CONSTRUÇÕES METALO-MECÂNICAS MAGUE, L.”” 
(S. E. MONIZ DA MAIA & VAZ GUEDES, L.”*) 
Alverca — Portugal 


a a. 
o mg 


—— = - = “E tee 


Guindaste eléctrico de 8t/10 m em funcionamento no cais da Companhia Cimento Tejo 


PROJECTO E FABRICO DE 


PONTES ROLANTES 
GUINDASTES 
GRUAS. DERRICKS E GUINCHOS 
APARELHOS DE ELEVAÇÃO ESPECIAIS 
CONSTRUÇÕES METÁLICAS PESADAS 


TECNICA XI 


para impulsos com 
tempos de crescimento 
rápidos e grande amplitude, 


impulso de disparo 


impulso prévio 


E impulso com tempo 
a ca 


simples ou duplos, com sm E ee 


potências de saída até 1,3 Watts. 


novo. 


Vendas e assistência técnica em todo o Mundo 
Philips Portuguesa S.A RL, Rua Joaquim Antonio d Aguiar 66 LISBOA 
Philips. Dept? EMA, EINDHOVEN, Holanda 


TECNICA XIV 


impuíiso duplo 
com tempo de crescimento rápido 


10V - 100V e 


impulso de grande amplitude 


impulso duplo de grande amplitude 


GERADOR DE IMPULSOS DUPLOS TIPO PP 1122 


Impulsos com tempos de crescimento rápidos: polaridade positiva 
ou negativa; amplitude de 10. mV até 10 V; impedância de saida 75 
ohms; tempo de crescimento < 20 ns; impulsos simples ou duplos. 


impulsos de grande amplitude: polaridade positiva; amplitude 10 V 
até 100 V; impedância de saida 75 ohms até 666 ohms, consoante a 
tensão de saida; tempo de crescimento < 100 ns; impulsos simples 
ou duplos. 


Largura do impulso: ajustável entre 50 ns e 50 ms. 


Frequência de repetição dos impulsos: internamente, de 10 c/s a 
| Mc/s; externamente, simples a 2,5 Mc/s. 


Impulso prévio: disponivel separadamente; largura do impulso 
100 ns; tempo de crescimento 20 ns; amplitude 1 V até 10 V; 
impedância de saida 75 ohms. 


Atraso: (do impulso principal ou do segundo impulso num impulso 
duplo” ajustável entre 100 ns e 10 ms. 


“= PHILIPS 


instrumentos electrônicos 
de medição 


Integração, atribuições, perfuração e inscrição. Es- 
tando programadas as diversas partes doprograma, 
têm que estas ser reunidas, perfuradas e inscritas 
na máquina. Distinguimos, para maior clareza, a 
leitura para a memória da máquina, que se faz 
por uma ordem de leitura do programa inscrito, 
da incrição inicial que deve ser feita por ordem 
externa. Esta ordem é dada, como se viu atrás, 
compondo no painel um endereço (E) e premindo 
o botão (L): uma palavra que esteja perfurada 
em fita à cabeça de leitura da máquina irá ser 
inscrita na memória que tem aquele endereço. 
Representá-la-emos por 


INSCRICAO E 


A inscrição, palavra a palavra, dos dados e do 
programa é inconveniente e pode ser substituída 
pela inscrição de um programa de leitura, seguida 
de execução desse programa. Suponhamos que o 
programa de leitura deva ocupar às memórias de 
PRIL (princípio de leitura) a FIML (fim de lei- 
tura). A ordem de leitura 


L) LEI PRIL FIML 


inscrita por 
INSCRICAO L 


promoverá a leitura do programa de leitura, perfu- 
rado e posto ao leitor, quando a máquina fizer a 
execução de L). A execução de uma ordem, sem 
ser em curso de programa, também se faz pela 
ordem exterior dada no painel que indicaremos 
com 


EXECUCAO L 


Neste caso o programa de leitura (onde come- 
çaremos a fazer as atribuições dos endereços 
efectivos) será 


11) LEI 301 305 
12) LEI 201 201 
13) LEI 101 110 


e a ordem de leitura daquele programa 


1) LEI 011 013 


com a ordem exterior de inscrição 


INSCRICAO 1 


Finalmente a fita perfurada terá 


LEI 011 013 
LEI 301 305 
LEI 201 201 
LEI 101 110 
PAR 
do; dj; dp; dy; as; 
X 
MUL 301 201 501 
SOM 501 302 501 
IMP | 501 501 
PAR 


e a inscrição será feita por 


INSCRICAO 1 
EXECUCAO 1 
EXECUCAO 11 


O cálculo, quando o quisermos efectuar, será 
comandada por 


EXECUCAÃO 101 


PROGRAMAÇÃO DO CÁLCULO DE POLI- 
NÓÔMIOS USANDO CICLOS. Trata-se de redu- 


zir o programa as instruções 


MUL 501 201 501 
SOM 501 (ai) 501 


repetidas as vezes necessárias. 

Uma rotina de cálculo aritmético é sempre 
constituída por uma sucessão de instruções a cum- 
prir. Queremos assim aludir ao procedimento do 
calculador humano que faz cada coisa (sempre 
bem definida) por sua vez, seguindo de um prin- 
cípio a um fim. O mesmo acontece com os com- 
putadores aritméticos. A diversidade de movimentos 
que um e outro fazem no decurso de um cálculo 
é, também, limitada. 

Por isto um programa pode ser representado 
pela relação rectilínea dos seus passos, e o próprio 
programa constituído pela lista das instruções 
que vai cumprir em sequência. As instruções em 
tais programas não podem, naturalmente, ser 
mais mumerosas que as memórias da máquina. 
Numa máquina de 1000 memórias, quando cada 
instrução levasse em média 1 milisegundo a 
cumprir-se, um programa rectilíneo duraria menos 
de 1 segundo. 
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A possibilidade de programas mais longos 
reside na existência de troços que se repetem no 
programa rectilineo, com modificações que podem 
automatizar-se fácilmente. Durante o cálculo a 
máquina, convenientemente programada, executa 
repetidamente os ciclos de operações, comandados 
por uma zona fixa da memória, que correspon- 
dem aos troços repetidos do programa rectilíneo. 
Um programa com ciclos programados represen- 
tar-se-á com o auxílio de anéis que simbolizam 
os ciclos. Um ciclo de instruções num programa 
é pois uma sucessão de instruções cujo cumpri- 
mento o computador repete. 

A programação de um ciclo ocupa-se geral- 
mente das seguintes partes: preparação, cálculo, 
decisão e modificação. Na determinação do valor do 
polinómio, o cálculo é a multiplicação por x se- 
guida da adição de aj, que representámos sim- 
bolicamente por 


M)  MUL P X P 
Sd) SOM PP A P 


onde P representa a memória de cálculo onde 
acabará por inscrever-se o valor do polinómio, 
X a memória onde está x e A; a memória onde 
se tem o coeficiente que intervém no momento. 
Cada vez que executa as operações de cálculo 
o computador tem de reconhecer se o valor do 
polinómio já está obtido (caso em que deve pas- 
sar a outra parte do programa), ou se deve repe- 
tir o cálculo. Chama-se decisão a esta parte do 
ciclo. Quando deva repetir-se o cálculo, algumas 
instruções deverão ser previamente modificadas ; 
é a modificação. A preparação de um ciclo é cons- 
tituída pelas instruções que sejam necessárias 
para pôr os dados na situação que se supõe no 
início do ciclo. 

Assim, neste caso, a decisão é programada 
com a instrução de contagem e a de comparação 
e transferência: 


CONT) SOM (1) 
COMP) SAL (4) 


CONTA CONTA 
CONTA FRENTE 


A modificação consistirá em substituir na ins- 
trução 


Si) SOM P Ai Pp 


o endereço A; por Aj;1. Será conseguida pela 
instrução 
MODE) 


SOM Si (1000) S; 
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Finalmente a preparação será constituída pelas 
instruções (se forem precisas) que colocam as 
em P, O em CONTA, as constantes 4 e 1000 em 
memórias convenientes. É evidente que estas 
inscrições se podem fazer todas por leitura. 
Todavia a passagem do conteúdo de A, para P 
põe um problema que aparece frequentemente. 
Uma instrução. 


SOM Ao (0) P 


coloca em P o conteúdo de A, com a dificul- 
dade de utilizar um zero que tem de estar ins- 
crito. Pode haver vantagem de pôr permanente- 
mente em zero o conteúdo de uma das 1000 
memórias; a memória 000, por exemplo. Ou, de 
uma maneira ou outra, criar uma instrução de 
limpeza 


LIM B C 


que anule os conteúdos das memórias B, 
B+-1,... C. Ou, ainda, de criar uma instrução de 
translação 


TRA A B C 


que transfere os conteúdos das memórias de A 
a B para C e seguintes. 

Vamos supor que se escolheu para o computa- 
dor ABC a solução que consiste em criar os 
circuitos para estas duas instruções a que se dão 
os códigos em linguagem máquina 


8 (a) (b) (c) 
9 sas (b) (c) 


Exercício 6. Escrever em linguagem máquina 
um programa compacto para o cálculo e perfu- 
ração do valor de um polinómio de 4.º grau, de 
que são dados inicialmente os coeficientes, para 
um valor do argumento, dado também inicial- 
mente. 


Exercício 7. Escrever em linguagem máquina 
um programa compacto para o cálculo e perfu- 
ração do valor de um polinómio de grau que 
não excede 10 para um valor de x que será 
dado, juntamente com os coeficientes, em tempo 
oportuno. 


SUMÁRIO 


Este artigo é uma versão simplificada, para fins pedagógicos, do que seria um curso de 
programação de um computador real. Aproveita-se a oportunidade para esclarecer algumas 
das mais importantes questões postas por quem se aproxima de um grande computador. 


SUMMARY 


This is an introductory paper about the main points of the big electronic computers, 
The author has two ains: to be as simple as possible and to answer as many of the funda- 
mental questions as possible. 
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NOTAS INFORMATIVAS 


C. D. U. 624.311.5/94 


Elementos sobre a produção e o consumo de energia 
na rede eléctrica nacional 


— Elementos extraídos das estatísticas mensais do Repartidor Nacional de Cargas (R. N. C.) 


Nora: As produções e os consumos das empresas do R, N. C. representam 
cerca de 94 */, dos totais do Pais, 


JANEIRO 


I — Breve nota mensal 


No aspecto hidrológico, o mês de Janeiro apre- 
sentou-se, no conjunto, muito abaixo da média. 


II — Elementos gerais (GWh) 


a) Mensais 
| 1962 | 1963 = | 

0 
Produção hidráulica (Ph) ...| 375,9| 404,5 | + 8 
Produção térmica (P+)..... 9,5 9,2 |+ 59 
Produção total (PT). ..... 381,7 | 413,7 + 8 
Exportações (Ex) .. ...... 16,8 1,2 )|— 93 
Importações (D) -. cc. e cus ss 0,0 0,0 0 
Produção para con- | 
sumos não perman. (Pcnp) 99,6] 79,0 + 33 


Produção para con- ' 
sumos permanentes (Pcp) (!) | 305,3 | 333,5, + 9,2 (3) 
Coeficiente de hidraulicidade | 159] 0,76 a 


NOTAS: 
(1) A produção para consumos permanentes (Pcp) é determinada 


pela seguinte expressão : Pcp = Pr — Penp 4-1 — Ex 


(2) O aumento percentual de produção para consumos permanen- 
tes, tendo em conta a incidência dos domingos e dias especiais, é 
respectivamente de 11,8"), e — 9!g. 


III — Diagramas de carga dos dias característicos 


| 4.º feira; 


| 1963 1964 

Produção hidráulica (Ph) MWh| 12503 14123 
Produção térmica (Pr) MWh 0 0 
Produção total (PT) MWh| 12503 14123 
Trocas com [ Export. (Ex) MWh U 0 
Espanha  Rmport. (1) MWh 0 0 


TOTAL 
Prod. para cons. perm. 


Pr-+(L-Ex) MWh| 12503 | 14123 
(Pcp) MWh| 10610 | 11474 


Prod. para cons. não perm, (Penp) MWh 1593 | 2649 

Potência máx. MW 692 761 
se Pr + (1-Es) Potência min. MW 292 348 
SETH OC | Utiliz.da ponta horas | 18,1 18,6 
ER Factor de carga 06 + 0,77 
5 | Potência max. MW 617 6143 
do p Potência min. MW 208 9291 
Õs ep Utiliz. da ponta horas 17,2 | 17,8 

Factor de carga 0,72 0,74 
TRONICA 


306 


TE TER 

II SIIISLIIS ISIS 

RRRRRESREREERHRES:. 
ER LAY 


IV — Energia armazenada nas principais albufeiras 


| No fim do mês 
Albufeiras : a 
GWh “o (3) 
Paradela. +. . cc. cc «| 221,4 99,9 
Venda Nova . + . «+» «| 124,9 97,6 
Salamonde . . «.. cc vv. 22,9 83,0 
Caniçada .ccccc e cr. 29,0 87,6 
Cabril i . . cc 00 0 o «| 2966 87,4 
Castelo do Bode. . ., ..| 148,6 88,1 
Guilhofrei . . .... . 6% 4,0 48,2 
Lagoa Comprida . +... .. 23,5 (2) 15,9 
Santa Lusa . vo.“ 4. 44,7 72,6 
Pracana .» co. cc... 9,5 73,6 


Póvoa os O W mid é O iss q 12,7 (3) 97,1 
Total. ..| 932,8 89,4 


Notas: 


(1) Coeficiente de enchimento em energia das albufeiras, 

(2) Inclui 1,8 GWh armazenados em Vale do Rossim no 
início do mês e 1,2 GWh no fim do mês, 

(º) Inclui 2,6 GWh armazenados no açude do Poio no 
início do mês e 2,6 GWh no fim do mês, 


e A TT e a e a e a o o e e SL 


PORTUGAL 


SIEMENS 


..... 
AMAM: 
40 nO" 
.* o o 


“e. 


MOTORES 


sea, . 
e... 


PORTUGAL 


SABUGO 


AG 


SIEMENS-SCHUCKERTWERKE 


FABRICADOS SOB LICENÇA DE 
BERLIN-ERLANGEN 
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SIEMENS COMPANHIA DE ELECTRICIDADE 
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TECNICA XV 


FÁBRICA NACIONAL DE FELTROS INDUSTRIAIS 
SOCIEDADE LIMITADA 
APART.: O ESCRITÓRIO E FÁBRICA 


TELEF.: 93 OVAR 
TELEG.: FELTROS PORTUGAL ESTRADA DE S. JOÃO 


FELTROS, FILTROS E TECIDOS INDUSTRIAIS PARA TODOS OS FINS: 


EM PEÇAS, PANOS, MANGAS OU SACOS, 
DE QUALQUER FIBRA 

PASTA PARA PAPEL TECELAGENS 
Bea =D 4Ê ACABAMENTOS 
FIBROGIMENTO TINTURARIAS 
CURTURMES LAVANDARIAS 
FIAÇÕES DE ALGODÃO ESTAGPARIAS 
410 NA 8 FILTROS 
INDÚSTRIAS QUÍMICAS DESPOEIRAGEM 


IMPERMEABILIZAÇÕES 


TORNA TUDO ESTANQUE 


Impermeabilização é vedação perfeita de: 


Terraços, Empenas, Algerozes, Clarabóias, 
Paredes húmidas, Coberturas, Fendas, 
Juntas de dilatação, etc., etc. 


Fácil aplicação — Elasticidade excepcional 
Não endurece nem apodrece 


R 1 ) F AT EX À base de borracha + alumínio 


Reforça a acção impermeabilizante do «Rubio», dando às construções um aspecto agra- 
dável — Protege todos os materiais contra a humidade, a ferrugem e a corrosão 


FLUIDO E MASTIC 


RUBIO E ROOFATEX dois produtos que se completam 


Representante exclusivo : (M A 


Rus Palmira, 46-A-B-C — LISBOA-1 — Telefs. 830163-8426 39 — End. toleg. «ORFER» 


TECNICA XVI 


C. D. U. 517.514.2: 536.7 


PROBLEMA DE EXTREMOS 


por ANTÓNIO GOUVÊA PORTELA 


(Conclusão) PROFESSOR DO |. S.T. 


Generalização ao espaço Rm 


0) O objectivo deste capítulo consiste em generalizar os tratamentos formais referidos ante- 
riormente a um espaço de n dimensões e ainda sugerir uma solução mecanográfica. 

O problema central reside no tratamento dos termos negativos de expoentes pares ou termos 
dispondo de alguns expoentes ímpares e que podem ser tornados negativos e determinar se esses 
termos podem ou não, em determinadas condições, ser maiores do que os restantes termos positivos 
pares ou termos de expoente impar que são positivos quando o termo ímpar é negativo. 

Resolvido este problema central para cada um dos termos negativos de expoente par ou ter- 
mos de expoente ímpar que podem ser tornados negativos, verificado que tal situação não ocorre 
então poderá afirmar-se que, no ponto considerado, a função tem um extremo (mínimo). 

O mesmo procedimento poderá empregar-se quando o extremo é um máximo, mas então 
deverá testar-se a negatividade da série. 

O método envolve tomar-se as derivadas até uma certa ordem; em física, interessa, quanto 
muito, estudar a série até à 6.2 ordem, pois não tem sentido físico, em geral, um extremo que se 
decide para além da 6.2 ordem. 

Formalmente o método é válido para qualquer ordem. 


1) Problema central 


ao! aon 

Sendo dados: um termo Mo > Mo hi... hy 
e 

ail ain 


m termos M; > M; hi... hy 


Determinar quais as condições em que: 


a0l aon 


Mo hi ... hn (1) 


sabendo que : 
Mo e M; são números finitos. 


ão; € ai são números inteiros finitos e positivos. 


2 
| = (1,2 .. in) 
j=(1,2...n) 
3) |h|<|9|] e |ô] se pode fazer tão pequeno quanto se desejar. 


Para encontrar as condições acima referidas, proceder-se-á à resolução de dois problemas. 


1.1) Estabelecer as condições necessárias e suficientes para que se verifique simultâneamente 
as m seguintes desigualdades : 


aúl aon 


4) | Mo hi .... hn 
e ainda 2) e 3). 


ail ain 


so M; h; e... hn 
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1.2) Provado que os sistemas 2), 3) e 4) são satisfeitos, propor em que condições a desigual- 
dade (1) também é satisfeita. 


2) Teorema I- À 


Para que o sistema de m inequações seguintes seja satisfeito, quando |h;]|<l|d| e [9] 
suficientemente pequeno, 


Mo hu.....h ” M; hu......ho (5) 
Será necessário e suficiente que: 
[9]. [xy] < [Ai] '6) 
para x; > |N| e |N| suficientemente grande. 
Sendo : 
«j = (ao; — aj) donde (2; ser um número inteiro) 
lh | = | q|* 


(7) o < |Ih|< |qgl<1 donde x>1 
"1 o Elalá ly 


o | 


| ainda 2) e 3) 


x ok 
Com efeito : 
Substituindo em 5) |h;| por |q|* teremos: 
2 ao) Xj 89) 26 
IM, | . Ig > IM. Igl j=1...... m 
e x; >|N| 
ou 
” Mi =À : 
é (ao; — ay) x; < En TR . (Ln |q|) j==1...cm 
0 
e x > |IN| 


Note-se que [gl << 1 donde Ln|q/|/<O 


ou substituindo por a; e A, téremos finalmente 


| o x <A; [2] |x] < [Ai] 
e ou (m, n) (n, 1) (m, 1) 
| y >I|NI e (8) 
x > |N| 


3) Discussão de desigualdade | [2] . [kj] << [Ai]. 
xy <|NI| 
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Para o efeito, examinemos todos os menore [2,1] da matriz rectangular [2;;], onde k 


representa p linhas e | p colunas (com p< [e] extraídas da matriz [2;]. 


Vamos designar por [/2».:] qualquer dos ( 4) menores distintos que podem ser formados 
conservando p linhas bem determinadas e combinando as j colunas p a p. 
Vamos admitir que (9) det [2x.:] = O para todos os ( É] menores distintos formados. 


Então, é porque existe uma ou mais linhas que são combinações lineares das restantes. 
Mas se det [fu:]=O0 para as p linhas consideradas, então existe um vecto [7] 


tal que: 
[7x] [2%,;] = [0] para as p linha referidas. (10) 


Sendo )x números finitos. 


Vamos impor ainda que yu FO k=1....p, isto é, só existe uma combinação linear. 
No casa de existir mais do que uma combinação linear, tratar-se-á separadamente cada com- 
binação. 


Seja q uma qualquer das linhas que entra na combinação linear; então: 
Vk sa , 
a pa . [ Xk, ] =— [ x, 1] para k — Pp (11) 
Ya q 


e designando por [2,;] o vector dos coeficientes da de inequação q ou seja da linha q da 
matriz [2]. 

Mas, dum modo geral, [2x,;]. [x] << [A«] e portanto também o será para as (p — 1) 
linhas que resultam das p inequações em estudo, suprimindo a linha qg. 


Donde 
< 
| - Lag] - [x] e ou Hg TA] (12) 


(1,p—?1) (p—1,n) (n,1) (Lp -1) (p—1,1) e k&q 


Quanto ao símbolo da desigualdade as seguintes hipóteses podem ser consideradas : 
a) e > 0 para todos os (p— 1) k. 
q 


Então a desigualdade (12) toma a forma: 


[| ta y<| E] [A] ek+aq (13) 
Ya Ya 


mas, por outro lado, será, para a inequação q: 


Leg) e Li) ] Soa (14) 


e conjugando (13), (14) e (11) obteremos: 


A > [9] [il >| de] -[4 k + q (15) 


ora a desigualdade (15) é verdadeira para qualquer q, e portanto para todas as inequações que 
fazem parte das p inequações onde se reconheceu existir o) menores [2w,1] cujos determinan- 


tes são nulos. 
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Como Aq e Es - [Ar |] são constituídos por elementos finitos e fazendo x! = 7j.xXo 
Ya 


Xo, Xj E O 
e 
<< co e xo se pode fazer tão grande quanto se desejar, poderá (15) escrever-se 


onde 


da > [a] [4)>— [É] tal ço (16) 
Xo q /q Xo 


e para xo > | N | e N suficientemente grande será 


| q 


| Xo 


+ 


| «TA <tl e|3| 


tão pequeno quanto se desejar. 
Donde resulta que | [zj]. [xy] | << |il (17) 
q 


que define um hiperplano que designaremos de hiperplano de degenerescência cuja expressão formal é 
[z;]. [xy] =0 (18) 
q 


haverá tantos hiperplanos quantas inequações interventoras na degenerescência ou seja p. 
A expressão (17) mostra que o módulo de [x j!. [7%] está limitado por |9 | o qual pode ser 


q 
tão pequeno quanto se desejar; bastava para tanto escolher | N | suficientemente grande. 
Da expressão (15) pode ainda concluir-se que 


o <a Ji) + A <| 2] | | Me | + 4a (15 a) 


b) adk, < O para todos os (p—1) yx. e para um certo q. 
Ya 


Então a desigualdade (12) toma a forma 


Fm à al. Dil > [A] Ad] (19) 


Ta 


e é facilmente deduzível que: 
AM > [2] Dol 


-[2| .tad>[ 5]. ta 


Neste caso [2%] não está limitado em módulo mas apenas em ordem. Com efeito, a 2.º ine- 
q 


(20) 


quação do sistema (20) veio, quanto muito, reduzir o espaço de existência de [2;] . [x;]. 
q 


Note-se que esta conclusão é verdadeira apenas para a inequação q; para as restantes veja-se 
a hipótese c) tratada a seguir. 


c) E. S O conforme o valor de k e q escolhidos. 
ta 
Vamos designar por BE q 2 , respectivamente as fracções que são => 0 ou <, 0. 
7a Ya 
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Assim E >o donde Es [ow] [x] << [ea [Ax'] 


Xq Ya Ya 
Vk” Vk” Vk” 
e «o dona e] msi. bai & Ea [Av 
Ya g ta 
el =1 >0 donde [x;] [x] < Aq 
Ya q 


Nesta hipótese haverá desigualdades num sentido e outras noutro, mas em conjunto o sistema 
não impõe condições da forma | [x] [x] | <|2| a nénhuma das linhas (inequações) do sistema. 


Teorema II A 


ão! aon > ail ain 
Se (21) | Mo hu .... In ii E coa i=1, m, 
ao! aon a ail ain 
então (22) Mo hu +... hn a É M; ha .... In 


excepto se na matriz [aj — aoi] = [2] existirem algumas linhas para as quais se verifique a igual- 
dade [7x] . [aj — ao;)] = [0] e os yx todos do mesmo sinal, porque então a desigualdade (21) pode 
não ser suficiente para permitir concluir (22). 
As letras indiciadas têm o mesmo sentido e significado já definido no Teorema I A. 
Com efeito : 


A desigualdade (21) pode escrever-se 


— [03] . [x] + Ai 
(23) 1>I|ql | 


[ag] Db] <A i=L.m 


Ora já vimos que | 
x > |N| e |N| suficientemente grande 


ouseja —l[x).I|Ixy)+ A >O 


Em geral, — [x]. [x] + A; não é limitado superiormente e poderá fazer-se : 
(24) — [9]. Iy)+tA=|P|>o 
e (25) |P|=hLng.(Ln|q/)! e q<I1 
e substituindo em (23) teremos : note-se que: 
Ln à (Ln | q/)-! Lnp <O 
1>|ql Ln|lgql<o 


Somando as m desigualdades (23) obteremos : 


dá me ei] o EO) rá 
m>s |qg| 
p= 
m  Infi(Lnlg|)! q m 
ou (26) 1> 35 lIq| =— 5 fi 
j=l M iai 
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m 
Se impusermos g<1 será 3 à<m onde li> o 3 ti 
ar imo 


Ora a desigualdade (22) pode tomar a forma : 


m — [0]. [x3] + Ai 
Los 3 "q (27) 
i= 


e esta desigualdade é sempre satisfeita se 1 q.e.d, 
Tratemos agora a excepção referida no Teorema II A, 
A demonstração foi possível porque admitimos que era, em geral, satisfeita desigualdade (24). 
Porém se entre algumas desigualdades do sistema se verificar que: [yu] [25] = [0] e )k 


éra [Ad] + 

-Ya 

+ Aq e portanto — [x]. [x;)] + Aq está limitado superiormente e daí pode a desigualdade 
q 


(24) não ser satisfeita. 


todos do mesmo sinal, já vimos em (15) e 15a) que 0<— [5]. [x] + Aq < 
q 


5) Estado do caso degenerado 


O caso degenerado verifica-se quando entre um certo número de inequações existe para as 
correspondentes linhas de Matriz [2;;], relações lineares entre os seus elementos. 
Sejam as inequações i=1,...p, e (p<n,m) aquela para as quais se verifica 


[yw] [mj]=[0] ek=1,...p 
e além disso y« ou (yu — 0) para todos os kk, já vimos em (17) que 


| Lay] by] |<]ó| (17) 


e por mais pequeno que seja || haverá sempre um | N | tal que para x, > |N | a inequação 
17) é sempre satisfeita. E tal é verdade para todos os q desde 1 a p. 


Também já vimos que a desigualdade (22) aplicada aos p termos se pode escrever: 
(Veja-se 27)). 


p — [aj] [xi] + Ak 
(28) ps [0] a 
$ a [ki] | Mk 
ou Le |] “IM 
Ora para x, > |N| e | N| suficientemente grande | — Lai) «Ix]l<|ô9|l e |º| tão 
pequeno quanto se desejar donde 
| — Lai) [x] a = 
(29) |lal à | € [8] 


e | | tão pequeno quanto se desejar e finalmente 


p — [99]. [| 
(30) 2 Iql & 


Mk 
Mo 


% 


S 
z da 


Mk | 
— oa 


Mo 


< 


e || tão pequeno quanto se desejar. 
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O problema fica reduzido a verificar a desigualdade numérica : 


1 E (31) 
k=1 Mo 
p 
6) Antes de prosseguir, convém demonstrar que 5 Mi é o mínimo de função (32) 
k==1 o 
p — [aj] [xo] + Ak P 
= 35 |ql| cida quando | — [e] « Dyl | + O e ainda 2 au = 0 
k=1 k k=1 o 
e Xi = 00 
Com efeito, 
p — [23] » [xy] M 
ç- 3 (gp RE, DE (33) 
k=1 Mo 
A condição necessária ou de estacionaridade de 7 é: = 0 j=1,..n 
Xj 
ou seja 
do p M [2] . [x] 
Del-s IqlY E um O (34) 
OXj k= 1 Mo 
para) == 1,...n 
mas já vimos em (29) que para x > | A | e | A | suficientemente grande será: 
| la | o). im) q | < |d | e || tão pequeno quanto se desejar e daqui 
| p M — [a] . [x] | p M 
|-23 : al “a -|- Da 3 | < |ã| (35) 
| k=1 * | Mo k=1 k o 
e |ô:| tão pequeno quanto se desejar 
p MK 
donde se TB o) SS PR 4 (36) 
k==] k o 
” p Mk 
e % | tão pequeno quanto se desejar ou no limite 53 a Me = 0 (37) 
k=1 o 


que estabelece uma relação a verificar em cada caso. 

Falta mostrar que se a condição necessária for verificada, também o será a condição suficiente 
ou de extremo. 

Com efeito : 

A condição suficiente propõe-se deste modo 


E > = 
A? q = [dot [E] [dy]>o (38) 
Ox dx. 
mas 
p — [21] . [x] 
(2 -a. SE lg) * 
929 ne k=1 k Mo 
dx dxj d x1 
p ae 
os Mk Lx) [xs] 
= > a. as | k (39) 
k=1 E k Mo + rg 
e novamente será : 
PO. — [08]. [x] p 
Mk Mk 
3 ma). k mis À) 00 1 Ó 
(40) 2. a | q | a “o Me ||) <t 
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e |)s| pode ser feito tão pequeno quanto se desejar, fazendo x; > |N| e IN] tão grande 
quanto o necessário. 


mas (41) ê add a | ss é] p | | a [2] 
(n,n) (n,p) (p, p) (pn) 


Sendo É ad 
Mo 


| uma matriz diagonal e tendo em vista (41), (40) e (38) poder-se-á con- 


cluir que : 


q q j Mk 
[dx]! | = | ; [doa = | [ae]? [ont | | E 


| [25] - [dx;] < | 9 |] (42) 


e | | poderá ser feito tão pequeno quanto se desejar, fazendo x; > | N| e |N| tão grande 
quanto o necessário. 

Portanto bastará provar que é maior que zero a expressão abaixo, para provar que 4º 9 >0 
e portanto que se trata de mínimo 


Mk 


o 


(43) [dx]". [a | 


designemos por 


| Toy] . fal 
k 


dy; = Las) « [dx] 


então (43) poderá escrever-se: 
| - [dy] 


Mk 
e para provar que (44) é positivo bastará provar que a matriz | | 


o 


(44) [dy]'. | x 


Mk 


| é positiva definida. 


o 


Mk 


o 


> 0 
q.e.d. 


Ora a matriz referida é diogonal e todos os seus elementos são positivos, 


7) Portanto o estudo das p desigualdades nas condições referidas no parágrafo 5 resulta 


k My à 
na comparação de 3 — com a unidade. 
k = o 


M 

45) Se sz = > 1 o sistema não tem um extremo no ponto em estudo 

k= o 

k | Mk a à 

46) Se 35 '-—— | << 1 então o sistema tem um extremo no ponto 

k=1 | Mo 
47) Se s | Mk — 7 então o termo M, e os termos Mk (k==1,—...p) têm uma 

k=1 | Mo soma algébrica nula 


Neste último caso se a função tiver um extremo este não se reconhece até à ordem em que 
foi feito o desenvolvimento em série de Taylor da função cujo extremo se deseja apurar. 


8) Tem interesse aplicar a formulação deduzida ao caso de | n = 2 
m = 2 
p=2 
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O NOVO TEODOLITO 


com calagem 
zenital automática 


O teodolito com os aperfei- 
çoamentos mais modernos, 
permitindo medições mais 
fáceis, rápidas e mais exactas 


Automatismo duma concepção surpreen- 
dentemente simples: Prisma com líquido 
sem partes mecânicas, sem desgaste, sem 
desarranjos, sem reparações. 


Detalhes no prospecto Th 154 E 


WILD 


Modelo T1-A C/imagem invertida 
Modelo T1-AE C/imagem direita 


Solicitem catálogo ou demonstração 
ao representante exclusivo: 


VALO PORTODAL, Lo 


LISBOA — PRAÇA DAS ÁGUAS LIVRES, 8 S/L 6 —Tel. 681127 


TECNICA 


/ 


XVH 


aços vazados 
ferros meehanite 


mecânica MÉDIA E PESADA 
construções metálicas 


COMPORTAS 

TURBINAS 

CONDUTAS FORÇADAS 
RESERVATÓRIOS SOB PRESSÃO 
PERMUTADORES DE CALOR 


Fabricação de uma serpenilna para fornalha 
* Principais características: tubo de aço 
Cr-Mo- 3' sch. 80 x» Diâmetro da serpen- 
tina 1600 mm +“ Peso aproximado 3500 kg 
* Desenvolvimento do tubo aproxim, 200 m 


DIVISÃO METALO-MECÂNICA 
COMPANHIA UNIÃO FABRIL 


VENIDA 2 DE JULHO, 170 — LISIQA 2 + TELEFONES 670421 E 672061 + TELEGRAMAS: METALOFABRIL 


TECNICA XVII 


Para que seja possível existir uma degenerescência é necessário que: 


[7x] [24,1] = [0], Veja-se expressão (10) e tt o 
7a 
e como [ [2] [x] < [Ai] i=1e2 


x > |N| 


o a 
=—<o0 e —<0 


an É EE 


e a desigualdade [z;;| [x;] | A;| pode escrever-se sob a forma: 


+ «xi— BXIC AS 
L—ax—bu<CA!S 


A condição de suficiência propõe-se (Veja-se 37) 


! M 
“| — | + em dim j=1 
o Mo 
Mi Ma 
| =] e 644 | === = j=2 
o Mo 
a a M 
náo —>-=— mm g e as duas equações degeneram num só g | HD lp|l—|=0 
al 222 o Mo 


ou simplificando g | Mi |+M;=0. 
A condição de degenerescência estuda-se da comparação com a unidade de 


é [Ma [M, [Ma] [MA IML > 
k=1 | Mo Mo Mo| | Mo | < 
> 
Tanto a expressão | Mi | + | Ms | = | Mo | como g| Mil + | Mz: | = 0 coincidem com as solu- 


ções estudadas anteriormente como seria de prever. 


9) Vamos proceder à generalização de alguns teoremas e lemas já deduzidos para funções de 
duas variáveis independentes. 


Teorema HI-A 


Designando por (i) qualquer combinação dos n índices is... in, pode provar-se que: 
n ão) pre n a(i)j 
48) [Mo Mhj [> ZM lh 
j=1 j= 
Se fôr: 
a()j > do) 


n 

2i;jzl 

n n j=i 

ad a(i)) > A ao j 

= j= 
hi CSIiGiIS? jm. 

a(ij)j | inteiros e 

a(ojj | positivos |h;| <ô ed tão pequeno quanto se quiser 
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Demonstração: 


A inequação 48) pode escrever-se: 
Mes | n a(i) j — aoj ) 
Mo El | hy | 


a(ijj— ão; > O 


oo 
49) Tx Ê 
(D=(0) 


Com efeito, a circunstância de R 
2 (ayj—aoj) > 0 
j=i 
a 
exclue a combinação a(ijj=ao; porque 2 (ao; — ao) =0, o que contraria a segunda desigualdade 
j=i 


n (6) 
invocada e o valor deste termo é 1] |h;| =1 
j=i 


Mai 
Mas seja Q = REA | para todo e qualquer combinação (i); Q é finito porque Ma É 2 
o 
e Mg =E 0. 
Ma 
Substituindo E por Q no segundo membro de desigualdade 49), obteremos: 
o 
co | a(ijj — ao j oo n a(ijj — ao) 
2 | Mei (tm -1)<o É H | hy | ue 7 
(D=(0) | Mo j=] (D=(0) A j=l 
mas 
co n a(i)j— ao j n oo a(i)jj— aoj 
3 DU (hil = 1 2 |hl 
(D=() j=i j=1 (i)=(0) 
e 
oa a(ijj — aào j É =| by) O 1 
> h; = — E ="— e h 1 
(i)j=0 hi 1—|h;| 1—|hjl | < 
Donde 
Q 2 (h my ) E | 
(ile=to) Mi=1 e = 1—|hyl 
epara |h;!< 3 e suficientemente pequeno será 
n 
—L get! 
=" 1—[|hy;l 
e |:| será suficientemente pequeno para que 
Qo Ji q.e.d. 
Teorema IV- A 
| o E a a(i)) 
(50) 2 Ar h > 2 |My| D- hs! 
k=l (1) | 
Se se verificar : 
| A;DX0 para k=1,...n 
p inteiro e positivo 
a(yj > 0 a A 
; inteiros 
(51) aWmk 2 2P 
n 
ZJam>2p 
j=1 
(o hi | << |9 e | | tão pequeno quanto se desejar 
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Demonstração : 


Do Teorema III A será; 


2p a(i)j 
(52) Ar he > 2 Mol À Ih] 
se for: aj > 0 
a > 2 e lh|< |9 
(53) e P | | | 
2 agi >2Pp 


j=: 


Se escrevermos as desigualdades do tipo (52) para k=1,2...mn tendo o cuidado de não 
repetir termos já incluídos nas desigualdades precedentes e somando ordenadamente teremos : 


2p co n a (i)) 


z Ar he > a | Me | À by] 
fu 


se se verificar (51) q.e.d. 


Nota: Daremos a designação de ordem de um termo da série ao número inteiro dado pelo 


somatório seguinte : n 
2 a()) 
j=1 
Teorema V-A 
n a(i)j 


Seja dada a série Sr = E My, Ih; 
= 0 =1 


Se 2p n  a(i)j n “p dg O 
(54) E Mi se h; e. Rod As hy e 3 ii RÉ RE 
ai a(ijj > 0 
> 


então Sr > 0 para |h;|<|9| e |7] suficientemente pequeno. 
(a(yy e p são inteiros). 


Demonstração : 


A série Sr poderá escrever-se: 


n a(i) j a(i) oo n a(t) 
(55) Sr= 2 Mod hy = > Mo, 7 h + ) My Uh 
()=o ()=o (1) j=1 
n = n 
2 am S2p 2 aiii = 2P 
Mas do teorema IV-A temos que: 
' a(i) " a(ji 
(56) Jam >S | My. Uh > mi Dhy 
se for 51) e como é 56), será: 
a(i) j ai)i 


O < TM lh > Mi h; 


(i)j=o (1) j=1 


e de 55) será finalmente Sr > 0. 
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Nota: Diremos que um extremo se decide na ordem 2p, quando 2p for a ordem da série 


finita e de menor ordem que verifica a desigualdade 54). 


Teorema VI- À 


A condição necessária, mas não suficiente, para que: 


2p n (id) n 2p 
: My Hh > 2 Ash (57) 
(1) j=1 k=1 
para lh/l| 0 ed suficientemente pequeno é a seguinte: 
nm a(iyk mínimo < 2p 
a(i) k minimo = 2 Pp e então será: 
58 ou Ka 
Em e it e 85) | a(i)k mínimo par (59) 
k 
o Mi >O 
sendo Ask > 0 e Pp inteiro positivo. 
Demonstração 
A desigualdade 57) se for verdadeira sê-lo-á também para h; = 0ej & k 
2P a(i) k 2p 
Donde É Ma, he > Ak hk para k = 1, n (60) 
Mas para satisfazer 60) duas hipóteses se podem verificar : 
Ou aí(i k inínimo == 2p 
” 2p a(i) k 2p 
e então 2 Mi he = Me mínimo < he 
(1) 
Donde Mi, mínimo >> Ak > O para satisfazer 60) 
ou ajyk mínimo < 2Pp 
a 2p à(i) k a(i)k mínimo 
e então à Ma hx o Mai) min, + he dd 
(1) 
Mai mín. a(i)k — a(i)k mínimo 
onde Sm] qu BS , hj 
o Mi 
e como anyk—awykmin > 0 
S= 146 
será | € e haverá sempre um [5] suficientemente pequeno, tal que |S| possa 
|he dl <ô6| 


ser mais pequeno do que um número arbitráriamente pequeno dado. 


Mas au;+ mínimo < 2p, donde se My min > O e aiy mínimo par, será satisfeito 60), mas se 


Mie) mínimo < O OU aç)k mia ímpar então 60) não é satisfeito. 
Portanto as condições 58 e 59 são necessárias mas não suficientes. 
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q.e.d. 


Teorema VII-A 


— Verificado que para todos osk =1,...n são satisfeitas as desigualdades referidas no teo- 
rema VII-A até à ordem 2 p. 

-— Verificado que para todos os termos de expoente par mas de coeficiente negativo ou para 
os termos de expoente ímpar (um pelo menos), as desigualdades do Teorema II-A se verificaram 
no sentido inverso, igualmente até à ordem 2 p. 

Então a série Sr é positiva definida. 


Demonstração 


Do Teorema I-A, temos que: 


n 


ar; m n 
61) [Me th ne 


aij 


|h;|< [3] e |3| suficientemente pequeno. 


Significando r qualquer dos s termos pares de coeficiente negativo ou ímpares (um expoente 
ímpar) e significando i==m os termos pares positivos da série até à ordem 2 p que existem na série. 
Tendo em vista a demonstração do Teorema II-A, a desigualdade 61) poderá ainda escrever-se : 


à O od * a aii y ' 
Me Mb <ç. 2 IM] e lh|<J|W|<l6] 
s n arj m n aii 
62) ou Z|iM |] |h| <Z|M|/Ulh| e lhi<l/9,| 
r=1 j=l i=1 j=i 
note-se que [|M |=M > 0 
e | aj= pares e positivos 
Do Teorema VI-A será ainda, por hipótese : 
3 aik gr 
(63) e Mi he > Arsh e A 20 
e Ih) <o 
e haverá k =-n desigualdades deste tipo. 
Donde m n aii nom aik n 2p 
(64) 2M 1H hj > Z 2 Mhy > 2 Arch: 
i==] j=l k==1 |=1 k==] 
e | bx | 0" 
mas (62) e (64) podem ainda escrever-se : 
m n i s n 
2. 2Ml0h > 2 IM: | HH lh 
Emi j=l f=1 j==] 
1 m n aij n 9 
o 2M0h > ZA 
2 jul i=1 ==] 
Donde m n aj 3 n arj n 2p 
(65) 2Mlh >> EM- WHl|h;l + 2 Akhk 
i=j j=1 s=i j=l k=1 
mas já vimos que : (Teorema IV A) 
sa 2p ce , a (i)j 
(66) 2 Arch > 2 [Mol .d [hi 
k=k () j=l 
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Substituindo (66) em 65 obtém-se : 


; n aii 8 n arj ou n 
2M0h > 2 IME lh] + & [Mol Elbal 
j= fee ==. 1 == 


a (i)) 


ou retirando os módulos será : 


m n aii 8 n arj oá n  a(i)) 
StT=5MHIh + EM Hh; + 2 My 1h > 0 
j=1 r=| ij=t (1) i= 


i=] j 


q.e.d. 
Conclusão 
É dada uma função v (Z1....Zn) e procura-se saber se tem um extremo no ponto 
P (Zw,... Zno), até à ordem 2 p e admitimos que 2 (Z1...Zn) é susceptível de desenvolvi- 


mento da série de Taylor na região do espaço que contém P no seu interior. 
Modo de proceder : 


1) Efectua-se esse desenvolvimento até à ordem 2p. 

2) Verifica-se a condição necessária referida no teorema VI A para todos k =1,...n. 

3) Verifica-se as condições do teorema VI A para todos os termos pares negativos ou ímpares 
existentes na série limitada até à ordem 2p. 

4) Então se 1,2 e 3) forem satisfeitos podemos afirmar que 9 tem um extremo em P, definível 


até à ordem 2p. 
5) Se as condições (nomeadamente 2 e 3) não forem satisfeitas poderá tentar-se uma série finita 


mais extensa de ordem > 2p e repetir a verificação. 


Finalmente tem interesse fazer uma referência ao emprego de ordenadores para resolver os 
sistemas de desigualdades lineares relativos ao teorema VI A. 
Com efeito, os métodos mecanográficos são particularmente úteis na resolução destes sistemas 


de desigualdades. 
Note-se que se as funções forem polinomiais, pode tratar-se desde início a série finita cuja 


ordem é igual ao grau de polinomio. 
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TECNICA XX 


C. D. U. 691,328:539.37 


REFLEXÕES SOBRE A FISSURAÇÃO 
E DEFORMAÇÃO NO BETÃO ARMADO 


LUÍS AROUCA 
ANTÔNIO MEXIA HEITOR 


Eng.ºS Civis (1.5. T.) 
2.4 Parte 


4. A INCIDÊNCIA DOS FENÔMENOS DE FISSURAÇÃO E DE DEFORMAÇÃO 
NO DIMENSIONAMENTO EM BETÃO ARMADO 


O estudo anteriormente apresentado, baseado na resistência do betão à tracção bem como 
na consideração dos efeitos de contracção e fluência, permite dimensionar peças de betão armado 
diminuindo apreciavelmente a incidência dos fenómenos de fissuração e de deformação. 

Para tal, há que ter em conta as reais características dos materiais empregados e complemen- 
tarmente, durante a execução da obra, é fundamental que sejam tomadas medidas adequadas para 
prever os fenómenos em causa. 

Referir-se-ão seguidamente alguns critérios de dimensionamento e disposições construtivas 
adequadas para o conseguir. 


4.1 Critérios de dimensionamento 


A consideração do fenómeno da fissuração no dimensionamento do betão armado pode ser 
encarada sob dois prismas : 


— dimensionamento de estruturas cujo condicionalismo exige completa ausência de fissuras, 
como é o caso de reservatórios para fluidos, estruturas existentes em meios agressivos, 
e outras. Nesta óptica, a multiplicidade de parâmetros envolvidos, como se infere da expressão 
(12), denota a irrelevância duma afirmação apriorística que inibe a utilização de aços de alta 
resistência de elevada tensão admissível ca, em que portanto a percentagem de armadura 7 
é pequena, em peças de betão armado isentas de fissuras. A garantia duma qualidade de 
betão que conduza a um elevado q'pr, parece mesmo, em princípio, assumir a posição de 
factor dominante. 

— dimensionamento de estruturas em que se admite um determinado grau de fissuração 
traduzido por uma largura máxima de fissura admissível. É o caso mais corrente de 
estruturas de betão armado, em que se despreza a contribuição do betão para a absorção 
dos esforços de tracção supondo as peças fissuradas abaixo da fibra neutra. Note-se; 
contudo, que mesmo que a largura máxima de fissura não ultrapasse os valores admissíveis, 
interessa limitar a extensão da gona fissurada da peça para impedir uma deformação 
excessiva. 


A análise do problema da fissuração que em si próprio tem um carácter estatístico, pode 
desdobrar-se na determinação de fórmulas empíricas para estimar o espaçamento das fissuras, 
e na determinação da relação entre o espaçamento das fissuras e a sua largura. 

Os parâmetros básicos que influem no fenómeno da fissuração do betão armado em função 
das condições de trabalho das peças são múltiplos e complexos, e o seu estudo tem sido abordado 
experimentalmente nos principais centros mundiais de investigação de engenharia civil, conduzindo 
à dedução de várias fórmulas, óbviamente de valor limitado, mas que, como diz Wastlund «servem 
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satisfatóriamente como prova, assegurando que a estrutura sé comportará bem no que respeita 
à formação de fissuras». Porém, este aspecto do problema sai fora do âmbito deste artigo. Interessa 
apenas focar a determinação da extensão da zona fissurada, independentemente da largura das 
fissuras, através da sua repercussão no fenómeno da deformação. 

Apresenta-se a seguir um exemplo em que é posto em relevo este fenómeno: 

Seja uma laje em consola com um balanço de 3,50 m numa pala de cobertura (elemento 
no qual o peso próprio assume uma importância principal). Convém, quer por condicionalismos de 
ordem estética, quer de ordem construtiva (escoamento de água e outros), que tal pala não apresente 
uma fissuração acentuada com a consequente repercussão na ordem de grandeza da flecha verificada, 
e como habitualmente, supõe-se que razões de ordem arquitectónica aconselham um dimensiona- 
mento com espessuras reduzidas. 

Considera-se um betão a que corresponde 7'br == 500 kg/cm”, que de acordo com as normas 
S. 1. A., conduz aos seguintes valores para as tensões de rotura em tracção por flexão 


após 7 dias 74=0,75>=3>=<y 500 = s50 kg/m? 
» 28 » íGp= 3>x< V 500 = 67 kg/em” 
» 90 » cp=1,15>x3>V500=77 kg/cm? 


Considera-se uma sobrecarga de 0,15 t./m*. Arbitrando uma altura no encastramento de 
15 cm e na extremidade de 6 cm, ter-se-á através da composição das diversas solicitações, um 
momento flector de lei parabólica, cujo valor no encastramento será 


3,5 3,5º 


2 3,5? 
Me = Em 0,15 + 2,4 | - u X< 0,06 + 


(0,15— 0,06) | = 2,02 tm/m 


Para que se não obtenham fissuras será necessário que o momento de fissuração esteja 
coberto em toda a extensão do balanço da pala. 
O momento de fissuração no encastramento é após 28 dias 


100 x 15º 
Mes = 67 >< a =— 2,51 t.m/m 


Note-se desde já que, como após 7 dias o momento de fissuração no encastramento é 
0,75 X< 2,51 = 1,88 t.m/m < 2,02 t.m/m 


não se deve descofrar e aplicar a sobrecarga. Já após 28 dias, o coeficiente de segurança à fissu- 
ração será de 


Se se entrar em conta com a fluência, então o valor de 9'p reduz-se (fig. 13) e por consequên- 
cia também o momento de fissuração. As normas S.1I. A. dão a relação 


c'p/0'br = 3/1,5 = 2 


e assim pela fórmula (5) vem £ ==5. Por outro lado, a resistência do betão em tracção por flexão 


dos 28 aos 90 dias varia com 7 == 1,4 (ver S.I.A. art. 63). 
Os valores de tensão de rotura em tracção pura são: 


pras = 1,5 V 500 = 33 kg'cm? 
cproo = 1,15 X< 1,5 V 500 = 38 kg/cm? 
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donde pela fórmula (28) tem-se aproximadamente 


33 + 38 5 9 
Hs (1,35 + 0,13 e =— 64 kg/em* 


em vez de 77 kg/cm?. 

Se a laje fosse descofrada após 7 dias sem se aplicar a sobrecarga, ela não teria fissuras visto 
que o momento de fissuração seria 1,88 t.m/m e o do peso próprio, 1,1 t.m/m. 

Porém, até à conclusão da fluência ter-se-ia pelas normas 5.1. A. (art. 63) o valor de 3=2,5, 
e o valor da resistência à tracção por flexão (após fluência) baixava para 61 kg/cm?. O momento 
de fissuração correspondente, 2,24 t.m/m, só daria um coeficiente de segurança de 


2,24 
2,02 


= 1,1 


Interessa agora comparar os valores que se obtêm para as flechas no caso da laje sem fissuras 
e no caso da laje com a zona de betão à tracção completamente fissurada. Os cálculos da laje não 
fissurada, fig. 14, (em fase 1) serão feitos com o valor do módulo de elasticidade do betão dado 
pelas normas 5.1. A. 


500 
500 + 200 


Ep = 600.000 = 430.000 kg/cm? 


55600 23500 12000 3550 

gomimo q 0413 0456 — ID40 2920 

metm/m) q 0875 1750 2625 350 

Fig. 14 — Laje em consola na fase | 
A expressão que dá a flecha é 
A M 1 AL 

is Mi q = | Mi. Mp. Tr] 
Eb I Ep 3 E 


em que M, é o momento devido a uma carga unitária aplicada na extremidade (secção em que se 
pretende conhecer a flecha) e Mp é o momento devido às solicitações. 

Para o cálculo do integral que aparece na expressão aplicou-se o método de Simpson, divi- 
dindo o vão da laje, 3,50 m, em quatro troços de 0,875 m. No quadro que se apresenta a seguir, 
é feito o seu cálculo 


ne | x ro fome [om | da 

O l do 600 O Ó Ó 
1 4 do 600 0,113 0,875 22 000 
2 a 23 500 0,456 1,750 31 600 
3 4 12 000 1,040 2,695 130 800 
4 1 3 550 2.090 3,500 24 600 
2 215 000 
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A partir das operações indicadas no mesmo quadro obtém-se o valor para a flecha de 


215000 >< 0,875 
3 > 430000 


É == 100 = 1,32 em 


Determinou-se seguidamente o valor da flecha para a laje fissurada (em fase II) partindo da 
hipótese que a armadura está adaptada ao diagrama de momentos flectores de tal modo que 7a é 
constante. 

A consola foi dividida do mesmo modo que no primeiro caso. Na primeira linha da fig. 15, 
indicam-se as alturas úteis h, na segunda está o valor h (1-Kx) e na terceira linha os raios de cur- 
vatura determinados de acordo com a equação (14) para 7a = 3000 kg/cm?. 

Na determinação das flechas pode proceder-se por etapes, entrando com os raios de curva- 
tura sucessivos no valor de 41º/2R, mantendo os pontos de tangência que descem de «.l em 


81 


virtude da rotação «== — .. 


Fig. 15 — Deformada da laje em consola na fase II 


Portanto f será dado pela sucessão 


A E q 3 1 E a s 
deita LR a a RA CR] 


/ 


Admite-se o último troço recto visto que a espessura não se reduz a zero e esse troço ficará 
no estado 1. 


A aplicação numérica dá 


a 3 2 1 1 
f == 100 X 0,875º | - jo =— 


imã | = 14,3 cm 
2,97 57 2,36 36 ZA 27 

Note-se, ante a diferença de valores obtidos num caso e noutro, respectivamente de 
1,32 cm para a laje consola não fissurada e de 14,3 cm para o caso de total fissuração, a extrema 
importância com que o fenómeno de fissuração se reflecte no da deformação. 

Claro que a amplitude da diferença registada se refere a um caso extremo, escolhido de pro- 
pósito para pôr em evidência o problema, mas de qualquer forma, significativo. 


Veja-se agora um caso mais corrente, que traduzindo em última análise a incidência da fis- 


suração na deformação, exprime igualmente a ponderação de toda a problemática analisada ao 
longo do artigo. 
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Seja, por exemplo, uma laje simplesmente apoiada com 15 cm de espessura e um vão de 
4,65 m (esbelteza 1/35) sujeita a uma carga total uniformemente distribuída de 0,735 t/m?. Consi- 
dera-se como tensão admissível no aço 7a = 4.000 kg/cm'. 
Resulta: 
Momento flector M = 0,735 X< 4,65"/8 == 1,98 tm/m 
Tensão no betão op = 110 kg/cm? 
Momento de fissuração Mr == 100 x 15*/115 = 1,95 tm/m 


Valores das flechas atendendo à fissuração e aos fenómenos de contracção e fluência do betão : 


Comprimento fissurado 
u = (4,65 — 0,56) / 4,65 = 0,88; Kx = 0,215 


Flecha estática 
Raio de curvatura R (expressão 14) 


R = 0,135 (1—0, e = 55,5 m 
Flecha (expressão 16) 
ERP: (1 — 0,88?) . 100 = 1,08 cm 
8 >< 55,5 


Flecha devida à contracção — és = 0,3 º% 
Raio de curvatura Rs (expressão 20) 


Rs = 0,135 (1 a — 00 490 m 
0,3 
Flecha (expressão 16) 
2 
fe=— —2º (40,88?) x 100 = 0,15 cm 
8 > 400 


Flecha devida à fluência — 9 = 2,5 


Determinação do raio de curvatura (expressões 21 a 27) 


= 1,98 — 0,292 
(100 < 13,5) 13,5 >< 0,93 > 400 100 
2 >< 100 
Pitt ao o da 


2 > 100 
De =1 +14 228100 10><0,292 — 9,30 


Ro =47 m 


Flecha (expressão 16) 
55,5 


Flecha total da laje parcialmente fendida 
f + fs + f; = 1,43 cm 

Flecha função do tempo 
És + fy == 0,35 cm 
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Note-se que o facto de se ter considerado uma peça com elevado grau de estabeleza e de se 
” “ pa . 4 == 9 ae 
ter dimensionado a peça com um aço de elevada tensão admissível (7a == 4.000 kg/cm”) não conduz 
a flechas cujos valores sejam superiores aos usuais no dimensionamento corrente. 


4.2 Disposições construtivas 


Conforme já referido, e complementarmente com um cuidadoso critério de dimensionamento 
que atenda à problemática da fissuração e da deformação em betão armado, impõe-se que a execução 
da obra seja sempre acompanhada dum conjunto de disposições construtivas adequadas à consecução 
de tal finalidade. 

As principais causas que levam ao aparecimento de fissuras nas obras de betão armado, apesar 
da sua problemática ter sido devidamente ponderada num dimensionamento criterioso, são funda- 
mentalmente a cedência das cofragens e a contracção. 


4,2.1. Cedência das cofragens durante ou logo após a betonagem 


Se o escoramento não suportar o peso do betão, pode haver movimentos no início da presa, 
quando o betão ainda tem pequena resistência à tracção, ocasionando fissuras, por vezes em locali- 
zações aparentemente inexplicáveis. 


4.2.2. Contracção do betão 


Na fig. 16, apresentam-se curvas de contraccção do betão em função do tempo para casos 
típicos de peças correntes em obras de betão armado. A curva 1 mostra a contracção de um ele- 
mento bidimensional (pequena espessura), como uma laje, a qual está sujeita aos raios solares, 
A curva 2 mostra a contracção de elementos volumosos, cujo endurecimento é dificultado por 
condições geométricas e em que a sua situação na estrutura os expõe menos à acção dos raios so- 
lares. À curva 3 mostra a variação no tempo das resistências do betão à tracção 7'br. 

Por exemplo, num caso corrente duma laje com os bordos livres, como seja uma consola ligada 
monoliticamente a uma viga de fachada, o elemento volumoso — a viga — tem uma exposição total- 
mente diferente da da laje, sendo a sua presa mais lenta. Como o aumento de resistência do betão 
à tracção não acompanha a variação diferencial de contracção dos dois elementos, há a possibilidade 
de fissuração. 


— 


E GE E ESET a 
Le obr Betão novo 


Betão velho 
Fig. 16 — Variação da contracção e da resistência Fig. 17 — Fissuração devido à diferença 
do betão à tracção com o tempo entre a velocidade de endurecimento e de contracção 


Destas diferentes condições de contracção, resulta, como se mostra por exemplo na fig. 17 a, 
que no caso de uma laje que foi betonada por cima de uma superfície de um material não sujeito 
à contracção mas com aderência ao betão, o encurtamento devido à contracção cs cresce com maior 
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rapidez do que a resistência do betão à tracção 7'br , ocasionando uma fissuração no betão novo. 
A fig. 17 b mostra o caso de se betonar com um betão novo (1) por cima de um betão que já está 
num estado adiantado da sua contracção (la). O desfasamento em tempo é igual a ti. Também agora 
a contracção :s atinge um valor importante numa ocasião em que o novo betão tem um baixo valor 
de resistência à tracção 7'br. 

Estes fenómenos não podendo ser totalmente evitados, podem ser atenuados se se mantiverem 
permanentemente húmidos, sobretudo o betão dos elementos de grande superfície, durante um pe- 
riodo de tempo razoável para que 7'pr possa aumentar com o tempo até que o betão tenha atingido 
uma elevada resistência à tracção. 

No caso das figs. 17 ae 17 b a situação ainda pode ser melhorada se molharmos intensiva- 
mente a superfície do betão já executado antes de se efectuar a nova betonagem. 


CONCLUSÕES 


Procurou mostrar-se ao longo deste artigo a influência dos fenómenos de fissuração e defor- 
mação no cálculo do betão armado e a forma como os seus efeitos devem ser considerados no dimen- 
sionamento em ligação também com os fenómenos da contracção e fluência do betão. Verificou-se 
ser possível dimensionar em betão armado com completa ausência de fissuras, o que, contudo, 
podendo conduzir a soluções menos económicas, só será de considerar em casos em que tal condi- 
cionamento seja exigido pelas condições de utilização das peças. Nos casos correntes, é perfeita- 
mente admissível um determinado grau de fissuração embora convindo que as zonas fissuradas não 
sejam muito extensas para não agravar excessivamente o fenómeno da deformação. 

Focou-se ainda o primordial aspecto de execução das estruturas que condiciona o fenómeno 
da fissuração, por tal forma que se afigura inútil qualquer dimensionamento teórico que não seja 
acompanhado de adequada execução em obra. 


NOMENCLATURA 


Na simbologia apresentada procurou seguir-se na medida do possível os critérios definidos 
pelo «Comité Europeu do Betão». 


b — largura genérica duma secção 
bo — largura da alma duma viga em T 
f — flecha genérica duma viga 
g — distâncias do apoio à primeira fissura numa peça solicitada à flexão 
h — altura útil duma secção trabalhando à flexão 
ho — espessura da aba duma viga em T 
ht — altura total duma secção trabalhando à flexão 
o N 
bh; 9'br 
E; : Ei a 
n = E, — coeficiente de homogeneização aço-betão 
q — carga uniformemente distribuída numa viga 
x — distância da fibra neutra à fibra mais comprimida duma secção traba- 
lhando à flexão 
y — distância da fibra neutra à primeira fibra mais traccionada, do diagrama 
de tracção no betão duma secção trabalhando à flexão (v. fig. 2, pág. 692) 
y — distância da linha neutra à fibra mais comprimida numa secção sujeita à 


contracção (v. fig. 10 pág. 701) 
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p—— 
D edi rs 
Do =1+ Va + - 
ny 
Ep — módulo de elasticidade do betão 
Em — área do diagrama de momentos solicitantes, na zona duma peça solicitada 
à flexão, trabalhando em fase II (v. fig. 9 pág. 700) 
Fe — área do diagrama de momentos efectivamente absorvidos pela armadura 
na mesma zona (v. fig. 9, pág. 700) 
I — momento de inércia duma secção solicitada à flexão 
x 
Kyx h 
E sá 
K, = 7 
L — vão duma peça solicitada à flexão 
Ma — momento flector diferencial de equilíbrio numa secção sujeita à contracção 
(v. fig. 10, pág. 701). 
Mr — momento de fissuração duma secção 
N — esforço axial (o esforço de tração tem sinal positivo) 
Np — força total absorvida pelo betão na zona comprimida de uma peça à flexão 
N'p — força total absorvida pelo betão na zona traccionada de uma peça à flexão 
R — raio de curvatura duma peça solicitada à flexão 
W — módulo de resistência duma secção solicitada à flexão 
Z — força exterior de tracção numa peça de betão armado 
Zr — força exterior de tracção que provoca o aparecimento de fissuras numa 
peça de betão armado 
É te — coeficiente de curvatura 
Eb 
— Aa — percentagem de armadura de uma secção 
Ab 
SER X 
io 
E 
Uy = h 
Ea — encurtamento devido à contracção de uma peça de betão armado 
Eb — extensão genérica no betão 
ti — extensão de rotura do betão em tracção 
Es — encurtamento devido à contracção duma peça de betão simples 
6 — ângulo de rotação do apoio duma peça solicitada à flexão 
s 
"La 
p= e 
ht 
b 
Mi o — em que b— largura da laje interessada (viga T) 
Ta — tensão genérica no aço 
Tao — tensão genérica no aço ao iniciar-se a fluência 
Sp — tensão genérica no betão 
Tb — tensão genérica no betão ao iniciar-se a fluência 
Cbr — tensão de rotura do betão em compressão (cubos) 
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